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NATIONAL  GAS  TURBINE  ESTABLISHMENT 
*An  introduction  to  the  engine  response  rate  problem 

PART  I 

Control  system  theory  relevant  to  the  engine  problem 


-  by  - 


H.  Taylor 


SUMMARY 

As  the  flight  conditions  of  gas  turbine  engined  aircraft  become  more 
severe,  engine  control  system  complexity  will  be  increased  and  the  perfor¬ 
mance  demanded  of  the  control  system  in  respect  of  the  number  of  parameters 
to  be  controlled,  the  range  of  control,  and  the  rate  of  response  and 
sensitivity,  will  become  more  difficult  to  achieve. 

An  increasing  effort  is  accordingly  being  devoted  to  control  system 
design  problems  and  to  consideration  of  the  power  plant  comprising  the 
engine  and  its  variable  geometry,  together  with  the  control  system,  as  a 
complete  and  integrated  unit. 

Since  the  engine  element  itself  is  a  dynamic  system,  the  response 
of  which  is  a  function  of  basic  engine  performance  and  operating  conditions, 
the  study  of  engine  response  behaviour  under  controlled  conditions  must 
result  in  the  theory  and  technique  of  control  system  and  engine  performance 
becoming  closely  related. 

In  Bart  I  of  this  paper  an  introduction  is  given  to  the  fundamentals 
of  control  system  theory  relevant  to  the  problem;  this  i3  intended  as  a 
guide  to  the  non-specialist  rather  than  as  a  treatise  on  servo-3ystems. 

The  theory  and  technique  of  servo-system  design  are  highly  developed 
and  provided  the  engine  problem  can  be  rationalised  by  acceptable  means, 
the  wealth  of  mathematical  technique  already  available  should  enable  the 
various  problems  to  be  resolved  by  established  methods. 

In  Bart  II  of  the  paper  (R.188)  the  engine  problem  is  considered  and 
an  analysis  made  of  current  literature,  much  of  it  irork  undertaken  in  the 
U.S.A.  by  the  N.A.C.A. 

It  has  been  established  that  the  response  of  the  engine  element  is 
basically  that  of  a  simple  exponential  lag  system  and  that  within  certain 
limits  the  engine  time  constant  (which  completely  defines  such  a  system) 
may  be  calculated  from  engine  component  performance  data.  The  effect  of 
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altitude  on  engine  time  constant  may  also  be  calculated  and  is  found  to  be 
considerable.  It  has  been  shown  however,  that  although  t'hiB  data  enables 
a  preliminary  assessment  to  be  made,  secondary  effects  including  additional  lag 
terms,  which  may  predominate  in  importance  in  relation  to  high  gain  fast  response 
systems  are  not  induded,  and  the  response  under  high  altitude  conditions,  when 
significant  changes  occur  in  component  performance,  has  yet  to  be  investigated. 

The  engine  experiments  required,  thus  relate  to  investigating  the  depar¬ 
ture  from  simple  theory  which  occur  in  the  engine  element  and  the  reasons  for 
them;  this  data  is  obviously  required  for  a  range  of  engines.  In  addition, 
data  is  required  which  will  enable  a  similar  assessment  to  be  made  over  as  wide 
a  range  of  altitude  and  forward  speeds  as  possible. 

Consideration  of  the  methods  for  calculating  engine  response  shows  a  need 
for  simplified  methods  which  can  be  applied  at  a  very  early  stage  in  future  engine 
development,  since  until  a  suitable  engine  test  facility  becomes  available,  the 
development  testing  of  control  systems  for  the  supersonic  high  altitude  case 
would  appear  to  be  very  difficult  to  simulate. 

Two  basic  methods  are  available  for  response  tests,  namely  step  change 
or  harmonic  input  functions,  and  it  has  been  shown  that  an  arbitrary  input  may 
be  substituted  for  the  step  changes  at  the  expense  of  computing  effort. 

Step  change  or  arbitrary  inputs  are  simple  to  apply  in  practice,  but  the 
information  obtained  is  mainly  limited  to  time  constant  and  gain  data;  further 
analysis  may  extend  these  results  to  yield  the  frequency  response  and  secondary 
effect  properties  of  the  system,  but  the  method  is  sensitive  to  experimental 
error.  For  the  development  of  high  gain  fast  response  systems,  there  would 
appear  to  be  no  alternative  to  the  more  complex  harmonic  response  method. 
Considerable  care  is  necessary  in  performing  these  tests  to  achieve  a  true 
harmonic  input  and  the  general  experimental  methods  used  are ' extended  to  their 
limit  in  respect  of  accuracy,  especially  in  relation  to  small  engines. 
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1*0  Introduction 


The  current  N.G.T.E.  programme  of  work,  which  is  considered  in 
Part  II,  is  concerned  primarily  with  the  response  of  the  basic  engine 
element  to  controlled  changes  in  fuel  flow,  but  this  must  be  undertaken 
with  full  consciousness  of  the  part  played  by  this  element  in  the  overall 
.  system;  the  form  which  the  programme  follows  is  dictated  by  the  use  to 
t  which  the  results  are  .  to  be  put  in  control  system  design.  As  knowledge 

is  gained  the  work  will  become  progressively  more  linked  with  servo-system 
design  and  with  the  overall  system  problems  of  sensitivity  and  stability. 

The  theory  and  development  of  control  devices  made  considerable 
strides  during  World  War  II,  but  for  security  reasons  the  work  was  not 
published  until  after  the  war.  A  wealth  of  data  and  specialised  technique 
was  released  eventually,  which  due  to  the  intensive  development  which  had 
occurred,  resulted  in  the  subject  tending  to  be  incomprehensible  except 
to  the  specialist. 

Since  the  engine  problem  is  only  partly  one  of  servo-systems, 
being  largely  a  performance  study,  it  embraces  a  range  of  interests  many 
of  which  are  unconnected  with  servo-system  design.  whilst  there  are 
excellent  books  and  papers  dealing  with  controls  and.  servo-systems,  there 
appears  to  be  a  need  for  an  elementary  paper  restricted  to  aspects  of  the 
subject  relevant  to  the  gas  turbine  response  rate  programme. 

Unfortunately  it  is  necessary  to  carry  this  further  than  justified 
by  the  present  state  of  the  response  rate  problem  due  to  the  style  of 
>  reporting  employed  in  current  reports.  In  the  case  of  many  American 

papers,  established  but  advanced  techniques  are  applied  to  problems  which 
are  in  fact  capable  of  presentation  in  a  more  elementary  and  as  far  as  the 
non-specialist  is  concerned,  more  intelligible  manner.  For  this  reason 
r  techniques  developed  for  control  system  synthesis  are  discussed,  even 

though  the  subject  in  general  is  outside  the  scope  of  this  paper. 

For  more  information  on  control  systems,  works  such  as  Theory  of 
Control  by  Macmillan,  and  Principles  of  Servo-mechanisms  by  Brown  and 
Campbell  (References  1  and  2)  may  be  consulted;  both  also  contain 
excellent  bibliographies  of  the  literature  relevant  to  control  systems, 
as  does  Reference  3. 

Harmonic  .response  has  been  introduced  in  this  paper  via  a  step 
function  approach,  since  the  latter  has  a  readily  understood  practical 
meaning  in  relation  to  gas  turbine  operation.  This  is  not  however 
intended  to  imply  that  the  step  function  approach  is  of  greater  signi¬ 
ficance  and  consideration  of  servo-system  design,  where  stability  is 
a  major  factor,  would  probably  favour  an  initial  approach  via  the  harmonic 
response  properties  of  the  system. 

Simple  control  systems  have  been  used  for  many  years  to  perform 
routine  tasks  and  to  control  remotely.  Controls  for  regulation  of 
*  engine  speed  and  industrial  processes  have  been  in  widespread  use  for  a 

long  time  and  over  a  number  of  years  the  theory  of  control  systems  has 
been  developed  independently  by  a  number  of  groups.  This  has  resulted 
in  common  problems  being  encountered  in  electrical  engineering,  comnuni- 
'  cations  engineering,  process  control  engineering  etc.  and  being  resolved 

independently;  inevitably  different  terminology  was  applied  in  the  dif¬ 
ferent  fields  and  interchange  of  technique  rendered  impossible.  As  a 
result  of  the  late  war  many  of  these  fields  were  brought  together  and  a 
common  approach  and  terminology  employed.  A  considerable  number  of 
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text  books  and  papers  on  the  subject  have  been  published  in  the  U.S.A.  and  these 
have  tended  to  place  emphasis  on  the  American  approach  and  techniques;  this  is 
especially  true  in  respect  of  the  Laplace  transform  method  which  is  universally 
applied  in  the  U.S.A.  to  simple  and  complex  problems  alike.  Whilst  it  enables 
discontinuous  functions  to  be  dealt  with  without  difficulty  and  is  a  powerful 
advantage  in  dealing  with  high  order  systems  and  for  synthesis,  it  is  not  greatly 
superior  to  the  classical  operator  method  for  the  analysis  of  simple  systems. 

Care  must  also,  be  exercised  in  the  solution  of  transient  problems  by  the  Laplace 
transform  method  to  ensure  that  the  "initial  oondition"  terms  are  correctly 
treated,  .  Nevertheless  much  of  the  published  literature  makes  exclusive  use  of 
this  method  and  it  must  be  considered  in  some  detail. 

Much  of  the  theory,  especially  that  related  to  harmonic  input  signals,  is 
closely  related  to  vibration  theory,  the  various  control  system  loops  being  analo¬ 
gous  to  the  elements  of  a  complex  vibratory  system,  the  stability  and  response 
behaviour  of  both  being  identical.  An  important  feature  of  control  system 
work  lies  in  the  importance  of  random  and  discontinuous  input  functions  and  in 
the  development  of  control  system  theory  as  a  tool  for  the  synthesis,  i.e.  devel¬ 
opment  and  improvement  of  particular  complex  systems,  as  opposed  to  mere  analysis 
of  system  behaviour'. 

2.0  Open  chain  and  closed  loop  control  systems 

Controls  divide  readily  into  open  chain  and  closed  loop  systems.  In  the 
open  chain  system  the  control  operates  to  provide  an  output  according  to  a  pre¬ 
determined  schedule  or  calibration  and  is  unaffected  by  the  resulting  output  or 
changes  in  output.  In  a  closed  loop  system  the  control  is  set  to  give  the 
desired  output  and  the  output  is  compared,  preferably  continuously,  with  the 
input.  In  this  way  control  action  is  taken  to  bring  the  output  to  the  desired 
value  and  to  counteract  any  tendency  to  disturb  it  from  this  value. 

A  rheostat  controlled  hot-plate  is  an  example  of  an  open  chain  system; 
the  rheostat  may  be  calibrated  in  terms  of  hot-plate  temperature  and  the  desired 
value  set  accordingly.  Any  change  in  calibration,  drop  in  voltage  etc.  cannot 
be  accounted  for,  i.e.  there  is  no  feedback  from  the  output  for  comparison  with 
the  command.  This  system  is  illustrated  as  a  block  diagram  in  Figure  1 . 

If  we  now  measure  the  hot-plate  temperature  and  compare  it  with  a  standard 
reference  temperature  controlled  by  the  command,  we  can  obtain  a  difference  sig¬ 
nal  of  appropriate  sign  indicating  whether  the  temperature  of  the  hot-plate  is 
high  or  low;  this  signal  may  then  be  amplified  and  made  suitable  for  operating 
a  motor  to  adjust  the  rheostat  and  so  bring  the  hot-plate  temperature  to  the 
desired  setting.  Any  change  in  oondition  will  result  in  the  generation  of  an 
error  signal  which  will  produce  a  corrective  action;  we  thus  have  a  closed 
loop  system  (Figure  1). 

3.0  Regulators  and  servo-systems 

The  terms  regulator  and  servo-system  often  lead  to  confusion  and  there 
is  probably  no  universally  agreed  definition  or  distinction  between  the  two. 

It  may  be  said  however  that  the  regulator  is  in  general  of  lower  performance 
and  is  mainly  concerned  with  maintaining  the  controlled  quantity  at  a  set 
level,  e.g.  the  conventional  voltage  regulator,  or  engine  top  speed  governor. 

The  range  of  control  is  usually  limited  in  respect  of  external  disturbances 
and  the  error  signal  may  generate  the  output  function  only;  this  will  usually 
result  in  a  steady  state  error  when  a  large  variation  in  input  quantity  is 
applied,  or  on  application  of  external  load.  The  latter  is  usually  known  as 
"load  droop. " 
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The  servo-system  however  will  normally  control  over  a  wide  range  of 
input  functions  and  will  control  with  random  input  or  load  variation.  The 
system  is  usually  power  amplifying  and  the  error  signal  will  normally  gener¬ 
ate  rate  of  change  of  output  in  addition  to  output,  thus  enabling  the  device 
to  exhibit  a  "follow  up"  characteristic  so  that  there  is  negligible  steady 
state  error  between  output  and  input.  This  is  probably  the  main  difference 
between  the  regulator  and  servo-system  but  must  not  be  taken  as  implying 
that  the  regulator  is  of  necessity  an  inferior  device,  since  in  many  appli¬ 
cations  it  gives  excellent  service  and  exhibits  a  high  order  of  accuracy. 

Finally  it  must  be  said  that  many  practical  systems  are  not  capable 
of  precise  identification  by  such  hard  and  fast  rules  as  the  above. 

There  is  less  ambiguity  in  the  definition  of  a  servo-mechanism,  which 
is  often  defined  as  a  servo-system  with  a  mainly  mechanical  output. 


4. 0  Basic  control  systems  and  their  behaviour 

4.1  Class  "0"  system  -  displacement  error  system  or  simple 
proportional  controller 


4.1.1  Open  chain  system 

Consider  a  positioning  device  comprising  a  long  shaft  free  from 
inertia  and  friction,  but  possessing  elasticity  in  such  a  way  that  it  has ' 
a  stiffness  k  in  torque  per  unit  angular  twist  units.  Also  consider  the 
shaft  subject  to  a  load  torque  Tp,  at  the  output  end  (Figure  2),  The 
symbols  used  are  identified  in  Appendix  II. 

Then  if  Q±  is  the  input  angle,  6±  is  the  angle  to  which  the 
output  would  set  in  an  ideal  system;  let  the  actual  output  angle  be  0o 


k  (6i  -  0O)  =  rL  • 

0i  -  0O  =  =  0  ...  ...  (l) 

where  0  is  the  error  angle  in  output  setting;  in  some  control 
system  work  the  opposite  sign  convention  is  used,  #iich  would  make 

—  negative  for  the  case  stated, 
k 

Thus  in  the  presence  of  a  load  torque  there  is  always  an  error  in 
the  output. 


4. 1 . 2  Closed  loop  system 

Consider  a  similar  system  to  the  above  but  with  the  addition  of 
the  feedback  loop  and  differential  device  providing  the  error  signal  0 
and  the  amplifier  of  amplification  factor  or  "gain"  p  which  is  necessi- 
*  tated  by  the  small  magnitude  of  the  error  signal.  0. 


k  (pe  -  0o )  =  rL 

By  definition:-  ®i  -  ®o  =  ® 

I'L 

H0i  -  p60  -  9o  =  Y 
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0o  a  .  6i  -  -L_  .  2l  . 

1  +  p  1  +  p  k 

e  =  _ 1 —  .  0i  +  ■— L—  .  . (2) 

1  +  |i  1  +  p  k 

Comparison  of  Equations  (l)  and  (2)  at  the  set  point  when  6-j_  =  o  indi- 

Pjj 

cates  that  with  the  closed  loop  system  the  error  has  been  reduced  from  —  for 
1  IV  ^ 

the  open  chain  to  - -  •  *r—  >  i* e*  for  an  amplifier  gain  of  10  it  is  reduced 

1  +  p  *■ 

to  one  eleventh  of  the  former  value. 

It  should  be  noted  that  in  the  absence  of  a  load  torque  the  closed  loop 

1 

system  exhibits  an  error  when  away  from  the  set  point  of  - -  .  0j_  ,  which 

would  class  this  system  as  a  regulator  rather  than  a  servo-meclianism. 

The  ircreased  lagging  error  with  increase  in  load  torque  is  an  example 
of  load  droop  and  is  a  common  characteristic  of  regulators. 

The  close  relation  to  feedback  amplifiers  will  be  noted  from  the  fore¬ 
going  and  it  is  not  surprising  that  much  of  the  knowledge  and  technique  employed 
in  control  system  work  was  developed  in  the  course  of  feedback  amplifier  design. 

4, 2  Class  "1 11  system  -  zero  displacement  error  system 

Whilst  consideration  of  the  Class  "0"  system  has  emphasised  a  number  of 
important  points  in  control  system  performance,  neglect  of  friction  and  storage 
terms  has  reduced  the  discussion  to  displacement  and  load  parameters  without 
the  conception  of  time  functions,  since  movement  of  the  input  produced  instan¬ 
taneous  movement  of  the  output. 

Consideration  is  now  given  to  the  closed  loop  system  discussed  previously, 
but  with  the  substitution  of  rate  of  change  of  output  proportional  to  error 
rather  than  output  proportional  to  error  (Figure  3) . 

The  conception  of  time  is  now  introduced  by  0j_(t),  0o(t) ,  0(t),  and 
rL(t)  as  time  functions. 

0(t)  =  0±(t)  -  0o(t) 

k  [k  J  0(t)dt  -  0o(t)]  =  rpi(t) 
k  (K  e(t)  .  iiaW]  .  law 

L  dt  -J  at 

,  x  d  00^)  ,  V  d  rL(t) 

k  K  0o(t)  +  k  =  k  K  0i(t) - T*p 


d  e0(t)i 

drL(t) 

dt  -J 

dt 

d  e0(t) 

= 

k  K  0i(t) 

d  rL(t) 

dt 

dt 

d  0(t) 

iC 

r  a  ©i(t) 

d  ri(t) 

dt 

dt 

dt 

Following  a  step  change  of  0^  ,  when  t-»ce  and  derivatives  are  zero:- 
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and 


Joss 


.0  =  0 


irrespective  of  the  value  of  0i  and  Tl  ,  where  0OSS  denotes  the  even¬ 
tual  steady  state  value  of  0O.  The  integral  term  has  thus  provided  a 
system  with  zero  steady  state  error  and  load  droop  for  a"  step  input  of 
displacement,  i.e.  a  servo-system  and  also  in  the  example  quoted  a  servo¬ 
mechanism. 


In  the  absence  of  load  torque :- 


K  0o(t)  + 


d  80(t) 


dt 


=  K  ©±(t) 


(3) 


Writing  the  operator 


D  s  t- 


d_ 

dt 


(D  +  K)  0o(t)  =  K  6i(t)  ... 

D  60(t)  =  K  e(t) 

If  a  step  input  of  0i  is  applied: - 

a  e0(t) 


(4) 


®i  -  ®o(t) 

®i  -  e0(t) 

When  t 

0o 

ec 

®o(t) 


-  k  at 


e  ?  e_Ki:  where  c  is  an 


arbitrary  constant. 


0 

0 


- .-»] 

angle/time 


1 


angle 


unit  of  time 


It  will  be  noted  that  K  has  the  dimensions 

If  we  write  t  =  y  ,  t  has  the  unit  of  time  and  is  known  as  the 
"time  constant",  to  which  considerable  significance  is  attached. 

Thus:- 

«o(t) 


04 


_t 
e  t 


which  represents  a  simple  exponential  lag  response  of  the  form  shown  in 
Figure  4. 


k 
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When 


When 


t  =  t 

&  ■  °-«2 
t  =  4v 
6n 

~  =  0.982 

8i 


* 


»■ 


which  gives  a  guide  to  the  response  time  for  virtual  completion  of  the  output 
motion  transient. 

Initially  when  a  step  change  of  01  i3  applied: - 

0  =  6j_  since  the  output  has  not  moved 

and  from  Equation  (3):-  d  0O 

Thus  the  tangent  to  the  output  response  -  time  curve  at  t  =  0  intersects 

the  —  a  1  line  at  t  =  ?  (Figure  4). 

0i 

t 

We  consider  now  the  case  of  a  step  change  in  input  velocity  i.  e.  in 

d  9j 

dt  ’  . 


Equation  (3)  may  be  differentiated  to  give:- 


- +  v - r — 

dt  dt51 

dt 

Writing 

d  0o(t) 
dt 

=  u0(t) 

and 

d  0i(t) 

dt 

=  wi("k) 

For  a  step  change  in  input  velocity  &>i  :  - 

,  t 

«i  -  «o(t) 


ece 


When 


.'.  From  Equation  (3) 


t  = 


Wn  = 


=  e° 


V 


and 


I 


Thus  <>io(t) 

e0(*) 

When  t 

0o 

•  •  Cj 

Hence  0o(t) 


In  the  steady  state  for  large  values  of  t:- 
®oss  =  ^i^  “  UfT 

The  response  to  this  step  input  of  velocity  is  shown  in  Figure  4, 
where  it  can  be  3een  that  after  the  initial  transient  condition  expressed 
_t 

by  the  tern  .  e  T  the  output  rotates  in  unison  with  the  input  with  a 

lag  error  of  oux.  This  may  be  minimised  by  keeping  t  to  a  low  value. 

The  term  Mix  is  often  known  by  the  misleading  term  "velocity  lag"  when 
in  fact  it  is  an  angular  displacement  lag  resulting  from  the  velocity  input. 
It  will  also  be  noted  that  when  the  input  is  arrested  the  output  "follows 
up"  and  cancels  the  lag  error  and  also  that  this  lag  is  present  without 
there  being  inherent  lag  in  the  controller,  i. e.  the  integration  has  had 
the  effect  of  adding  a  storage  unit  to  the  system. 

It  is  clear  from  a  study  of  the  Class  "0"  and  Class  "1 "  systems 
that  integration  within  the  loop  is  necessary  to  remove  steady  state  error 
and  it  is  of  interest  to  consider  a  double  integration  within  the  controller. 

4.3  Class  "2"  system  -  zero  velocity  error  system 

This  is  shown  in  Figure  3. 
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/  -t  \ 

=  m(l  -eT) 

=  wj,^t  +  're’rj+c1 

=  0 
=  0 
=  -<i>iT 

( 

=  Ujt  -  l  1  -  e  I 


k  K//e(t)  dt  .  dt  -  e0(t)  =  I^Ct) 


a?80(t)l  da  rL(t) 


k  K  8o(t)  +  k 


a8e0<t) 


dt* 


dt8 


=  k  K  9i(t) 


?  It(t) 

dt* 


With  aero  load  torque 


(D8  +  K)  60(t)  =  K  8i(t)  ...  . (5) 


T* 
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Hence  60(t)  s  ©i  (1  “  oos/K  t)  . .  (6) 

for  a  step  input  of  ©i. 

This  represents  a  continuous  oscillation  and  the  system  thus  requires 
some  form  of  damping  for  use  in  a  practical  case;  it  may  be  noted  in  passing 
that  the  simple  pendulum  is  an  example  of  a  Class  "2"  system. 

4.4  Order  of  control 


It  can  be  seen  from  the  foregoing  that  in  the  absence  of  time  lags, 
friction  etc.  controls  may  be  grouped  according  to  the  number  of  integrations 
within  the  loop  and  will  exhibit  distinctive  differential  equations  relating 
input,  output  and  error;  these  are  given  in  Table  I  for  displacement-displacement 
systems  and  it  can  be  seen  that  increase  in  the  order  of  control  tends  to  improve 
the  steady  state  behaviour  but  tends  to  diminish  the  system  stability.  It  is 
also  clear  that  satisfactory  operation  of  a  control  system  necessitates  a  care¬ 
ful  theoretical  assessment  of  the ' arrangement  if  the  accuracy,  rate  of  response 
and  inherent  stability  are  to  be  acceptable.  It  should  be  noted  that  there  is 
a  clear  distinction  between  order  of  control  and  the  order  of  the  differential 
equation,  since  these  are  not  necessarily  synonymous.  This  will  be  considered 
further  in  Section  4. 5. 1 . 

4.5  Consideration  of  a  practical  system 


The  simple  systems  discussed  previously  are  unlikely  to  occur  in  practice 
when  time  lags,  friction  and  inertia  effects  are  present.  These  systems  would 
in  any  event  be  little  use  in  practice  since  there  is  no  power  amplification  * 

and  the  input  and  output  cannot  be  physically  separated. 

4.5.1  Simple  system  with  and  without  time  lag 

* 

The  positioning  system  depicted  in  Figure  5  illustrates  a  basic  Ward 
Leonard  set  employed  as  a  position  control.  Friction  and  inertia  are  neglected, 
but  otherwise  the  system  has  the  basic  requirement  of  power  amplification  and  the 
ability  to  have  the  input  and  output  located  remotely,  A  rotary  potentiometer 
connected  to  the  output  shaft  provides  a  voltage  proportional  to  output  angle 
which  is  compared  with  the  voltage  provided  from  the  command  rotary  potentio¬ 
meter  measuring  input  angle:  the  resulting  error  signal  is  fed  to  an  ampli¬ 
fier  whose  output  controls  the  field  current  of  the  power  generator,  which  is 
in  turn  driven  by  a  constant  speed  motor.  The  servo-motor  is  fed  with  a 
constant  field  current  and  thus  has  an  input  voltage  proportional  to  error. 


Hence  the  voltage  applied  to  the  servo-motor  oc  6(t) 


Speed  of  the  servo-motor  oc  input  voltage  «  0(t) 


a  e0(t) 

dt 


(D  +  K)  e0(t) 
d  e0(t) 


k  e(t) 

K  [  0i(t)  -  0o(t)] 
K  ©i(t) 

K  0(t) 


* 


Hence 


-  15  - 


Report  No,  R.187 


This  may  be  compared  with  Equation  (4)  and  represents  a  first  order" 

system. 

Consider  now  the  system  depicted  in'  Figure  5  in  which  the  inductance 
and  resistance  of  the  generator  field  are  introduced,  i.e.  the  lag  of  the 
field  winding; - 

di.?(t) 

e  (t)  =  L— ^  +  R  if  (t)  =  |i*(t) 


if  (t)  =  0(t) 

i  +  t  D 


1  •  •••  •••  v7) 


Where  K,  =  — 


We  thus  have  a  simple  exponential  lag  interposed  between  the  error 
and  the  output. 


Hence  D  0o(t)  =  -  0(t) 

W  1+iB  U 


and  0o(t)  = 


Kq  0j(t) 

(t  if  +  D  +  Ko) 


•  •  •  • • •  • • •  (8) 


Where  =  K'  Kx 

and  K'  n  0(t)  =  D  (LD  +  R)  0o(t) 

For  a  step  input  displacement  0^;- 


0n  =  0-i 


0=0 

For  a  step  input  velocity  wj_;- 


0aa  =  -=■ 


Thus  the  addition  of  the  time  lag  has  not  altered  the  basic  first 
order  response  characteristic  of  zero  (displacement  error  for  step  input 
of  displacement  and  displacement  lag  -jg-  for  a  step  input  of  velocity. 

The  response  however  is  no  longer  a  simple  exponential  for  the  step  input 
case,  and  has  oscillatory  tendencies  governed  by  the  solution  of  the 
characteristic  equation  (tD8  +  D  +  Ko)  0o(t)  =  0  and  the  value  of  t 
and  Kg  (Figure  o).  Thi3  may  be  conpared.  with  Figure  4  for  the  basic 
system. 

Reference  to  Equation  (8)  shows  that  the  response  is  defined  by:~ 


D  (1  +  T  D) 
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The  governing  equation  with  a  second  simple  exponential  delay  would  take 
the  form:- 


0o(t) 


_ Jja _ 

D  (1  +  -ti  D)  (1  +  -ra  D) 


0(t) 


It  will  he  seen  later  that  where  the  lag  involves  an  inertia  and  spring 
force,  or  inductance  and  capacitance,  the  governing  equation  would  have  the  form:- 


"  D  (1  +  2£  ~  D  +  t4  DS) 

which  is  described  as  a  first  order  system  with  a  complex  exponential  delay,  the 
latter  tern  following  from  the  nature  of  the  roots  of  the  characteristic  equation. 


In  general,  lag  terms  have  a  de-stabilising  effect  and  it  will  be  noted 
that  the  order  of  the  differential  equation  is  of  course  raised  by  them. 


It  is  appropriate  at  this  juncture  to  consider  also  the  finite  delay  or 
time  lag  caused,  for  example,  by  lost  motion  in  a  system,  or  by  the  reaction  time 
of  a  human  operator. 


In  this  case  the  equation  for  the  delay  "element"  is:- 
0o(t)  =  0i(t  -  T) 

where  T  denotes  the  time  delay. 


Or 


e0(t)  = 


1  .  td  +  .  m 1 

2!  y. 


8i(t) 


-TD 


0i(t) 


■*.  6o(t)  =  ~ 


0i(t) 


TD 


/ 


The  finite  delay  is  thus  equivalent  to  an  exponential  operator  and 
introduces  in  effect  an  infinite  order  differential  equation;  it  follows  that 
the  finite  delay  may  have  a  severe  de-stabilising  effect.  This  is  a  most 
important  factor  in  a  practical  design  where  such  delays  are  frequently  present. 
The  effect  may  be  more  readily  appreciated  in  relation  to  the  harmonic  response 
locus  and  is  considered  again  in  1 2. 0. 


if,  5. 2  System  with  viscous  friction  and  inertia 

V 

The  systems  discussed  previously  are  unlikely  to  occur  in  practice  with¬ 
out  the  presence  of  friction  and  inertia  effects.  A  more  practical  remote 
position  control  (r.p.o. )  will  now  be  considered,  similar  to  the  last  case  in 
that  the  input  and  output  may  be  spaced  as  desired  and  considerable  output  power  » 
controlled  by  a  negligible  input  power.  This  system  is  depicted  in  Figure  5 
and  will  be  seen  to  comprise  a  constant  armature  current  motor  driving  the  output 
to  which  a  rotary  potentiometer  is  fitted  measuring  output  angle;  the  voltage 
representing  this  angle  is  compared  with  the  voltage  equivalent  of  the  input 
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angle  (as  set  by  the  input  potentiometer)  and  the  resulting  error  signal 
fed  to  an  amplifier  driving  the  split  field  winding  of  the  driving  motor. 

We  assume  that : - 

Motor  torque  rjj(t)  oc  field  current  if(t)  oc  6(t) 

d  ©0(t) 

For  viscous  friction,  friction  torque  oc  — — - 

Torque  absorbed  =  torque  applied 
d30o(t)  d  60(t) 

j  —■ Sl-i  +  f  — *  rL(t)  =  K0(t) 

dt  at 

Where  f  is  the  viscous  friction  coefficient 


e0(t)  = 


(J  D3  +  f  D)  6j(t)  +  fL(t) 
( J  D8  +  f  D  +  K) 

K  6j(t)  -  PL(t) 

(  J  D*  +  f  D  +'  K) 


k e(t)  =  d  (j  d  +  f)  e0(t)  +  rL(t)  ...  (1 

By  putting  D  and  powers  of  D  equal  to  zero  in  Equation  (9)  for 
the  steady  state;- 


0ss  *  K 


for  a  step  input  0^  and  constant  load  torque 

It  should  be  noted  that  the  system  is  a  servo-mechanism  in  that 
there  i3  no  error  due  to  input  disturbance;  the  effect  of  the  output  load 
to  produce  droop  is  a  characteristic  of  all  systems  in  which  control  torque 
is  proportional  to  error  only. 

For  a  step  input  of  velocity  and  constant  load  torque 


If  J  =  0  we  obtain  the  simple  first  order  system,  and  if  we  put 
f  =  0  the  basic  second  order  system  which  results  in  continuous  oscillation 
of  the  output.  Clearly  at  intermediate  values  of  f  and  dependent  upon 
the  values  of  J  and  K  there  are  output  characteristics  with  varying 
oscillatory  tendencies,  which  may  be  poorly  damped  but  are  not  of  course 
divergent.  When  f  =  0  the  system  oscillates  continously  at  a  mathe¬ 
matically  bounded  finite  value  and  is  said  to  be  neutrally  stable. 

.deamination  of  Equation  (11)  and  comparison  with  Equation  (8)  shows 
that  the  system  is  basically  a  first  order  control  with  a  simple  lag  added. 
The  differential  equation  however  is  of  the  second  order. 

The  exact  nature  of  the  solution  can  be  found  by  solution  of  the 
differential  Equations  (9)  and  (10). 
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The  solution  of  a  differential  equation  comprises  the  particular  integral 
representing  the  forced  or  sustained  solution  added  to  the  complementary  function 
representing  the  transient  or  unforced  system  response.  In  the  classical  method 
of  solution  a  forcing  function  must  be  established  and  the  particular  integral 
determined;  with  step  inputs  the  forced  or  sustained  solution  must  be  the  steady 
state  response. 

Hence : - 

Complete  solution  =  complementary  function  +  steady  state  solution 

The  complementary  function  is  obtained  by  setting  the  forcing  function 
at  zero,  giving  as  the  characteristic  equation  for  Equation  (9)  or  ( 1 0) :  — 

(J  Ds  +  f  D  +  K)  6(t)  =  0 

The  substitution  of  6(t)  =  em^  in  this  equation  gives 
emt  (J  m2  +  f  m  +  K)  =  0 

Hence  is  a  solution  provided  m  is  a  root  of  the  auxiliary 

equation: - 

Jma+fm  +  K  =  0 


mi,  ma 


Thus  in  the  absence  of  load  torque :- 


9(t)  =  Cj  e™!*  +  Ca  e®8*  +  0aa 


When  Ci  and  are  arbitrary  constants. 


When  f  =  0 


mi ,  ma 


- 


Now  the  basic  equation  form  for  a  second  order  system  was  given  in 
Equation  (6)  as:- 

0(t)  =  6i  cos/k'  t 

Hence  with  the  present  system  constants  in  the  absence  of  friction; - 

6(t)  =  0i  oo./}  t  which  represents  a  free 

oscillation  of  natural  frequency/^ 


system. 


Hence  /  —  =  %  where  «n  is  the  undamped  natural  frequency  of  the 


f  =  0 


mi,  ma  =  % 
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e(t)  = 

which  has  been  shown  to  be  oscillatory. 
f3  =  4JK 

,  ms  - 

and  0(t)  = 

which  represents  an  exponential  response  just  free  from  oscillation, 
hence  the  damping  is  said  to  be  critical. 

Thus  f0 

damping. 

We  can  now  write :  - 

5 

and  2£  % 

We  may  now  rewrite  Equation  (9)  In  the  absence  of  load  torque ; - 

(D3  +  2%  ofo  D)  6j(t) 

W  =  (D3  +  2£  (jfj  D  +  (On3) 

The  characteristic  equation  is  now:- 
(D3  +  2?  un  D  +  un3)  0(t)  =  0 

nb.,1%  =  -  gunljun/l-^ 

Where  %  / 1-2*  is  the  natural  frequency  of  damped  oscillation 
of  the  system. 

There  are  thus  three  forms  of  the  roots  of  the  characteristic 
equation  depending  upon  the  value  of  •  £: - 

X  <  1  Roots  conjugate  complex  with  negative  real  parts. 

System  underdaraped. 

X  -  1  Roots  negative  real  equal. 

System  critically  damped. 

X  >  1  Roots  negative  real  and  unequal. 

System  overdamped. 

Main  interest  is  oentred  in  the  underdamped  case  when; - 

r*i,  n*  =  -£  %  toi/1  - 
a  a  ±  Jb 


=  2 /JK  where  fc  denotes  critical 


actual  damping  _  f 


=  §  ,  since  wn  =  /§ 


0l  e+Jwnt  +  (^g-jwnt  +  0BS 

f_ 

“  2J 

(Ca  +  CUt)  eat  +  0ss 
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Hence  8(t)  =  eat  [cbe^  +  C4e-;3bt^  +  ess 

But  e#*  =  cos  bt  +  j  sin  bt 

Thus  0(t)  =  eat  £  (C!»  +  C4)  oos  bt  +  o  (Ca  -  C4)  sin  btj  +  0B8 

Since  the  solution  must  be  reed,  we  may  write 

0(t)  =  e”^^  ^  C®  cos  /(I  -  £a)  t  +  Qg  sin  Un  /<T  -  £a )  tj  +  eas 
where  C*  =  (Cs  +  C4) 

C6  =  5  (Ca  -  C4) 

The  constants  must  be  determined  on  the  basis  of  the  initial  conditions. 

No  purpose  is  served  by  considering  the  ensuing  algebra  here  and  the 
results  for  various  disturbing  functions  are  given  in  Table  II. 

A  typical  non-dimensional  plot  of  output-input  ratio  against  time 
for  various  values  of  £  is  given  in  Figure  6. 

Vilhen  £  un  t  =  4,  =  0.018  and  the  system  will  have  virtually 

reached  the  steady  state. 


A, 5. 3  Consideration  of  a  numerical  example 

Assume  J  a  49.10~*  slug  ft8 

f  =  0.157  ft  lb/rad/sec 

K  =  I.96  ft  Ib/rad 

The  value  of  J  would  correspond  to  a  motor  of  about  0.25  h.p. 
If  we  assume  a  moderate  input  velocity  of:- 
=»  15  rad/sec 

Then  Un  =  J~  =  /  —  =  20  rad/sec 

/  J  J  49.10  « 


and 


Time  to  reach  steady  state ,  t  =  = 

£«n 


0.25  sec 
1.2  rad 


0Ss  *  69° 


This  error  is  of  course  excessive,  and  la  due  to  the  system  being 
insufficiently  "stiff",  as  will  be  seen  later  when  considering  the  harmonic 
response  properties  of  this  type  of  system. 


Suppose  we  require  0a8  =  0.4  rad 
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=  0.4  = 


Hence 


=  5. 9  ft  lb/rad 


Wn  =  34.8  rad/sec 
£  is  now  rather  low. 


f 

vT 


For  a  fixed  value  of  J 
and  u-t 


Hence  if  we  double  f  and  K,  0SS  remains  unchanged,  but  £  is 
increased  by /2~. 


New  values 


49. 10~4  slug  ft8  ■  ■ 

=  0. 3"!  4  ft  lb/rad/sec 


ft  lb/rad 


=  0,4  as  before 


Ujj  =  49.0  rad/sec 

The  desired  performance  has  been  obtained  by  stiffening  the  system 
and  increasing  the  friction  coefficient,  and  it  is  instinctive  to  consider 
the  power  dissipated  in  viscous  friction. 


Power  dissipated  = 


f  w±a  746 


watts 


=  96  watts 

Since  the  motor  selected  was  only  rated  at  0.25  h.p.  i.e.  186.5  watts, 
more  than  50  per  cent  of  the  available  power  is  wasted  in  viscous  friction 
when  this  is  set  to  satisfy  the  specification. 

Thus  a  system  employing  viscous  friction  alone  for  damping  is  only 
of  use  for  low  power  ratings  and  is  inefficient;  there  is  also  the  danger 
of  large  changes  in  damping  occurring  due  to  changes  in  ambient  temperature. 

We  may  sunmarise  by  stating  that  in  the  practical  zero  displacement 
error  system  when  viscous  friction  and  inertia  are  present,  the  system  may 
have  three  forms  of  response  depending  upon  the  magnitude  of  the  system 
parameters,  namely  underdamped  when  £  <  1,  resulting  in  an  oscillatory  but 
stable  response,  critically  damped  when  £  =  1  giving  a  response  Just  free 
from  oscillation,  and  overdamped  when  £  >  1  leading  to  a  sluggish  non- 
oscillatory  response. 

The  achievement  of  a  suitably  damped  but  responsive  system  employing 
viscous  friction  alone  as  damping  leads  to  an  inefficient  system. 
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4, 6  The  achievement  of  a  damping  term  by  means  other  than 
viscous  friction  damping 

4. 6. 1  Output  rate  feedback 

Returning  to  Equation  (10),  in  the  absence  of  load  torque:- 

.  m  -  K  9l(t) 

0  [JD*  ♦  CD  «■  K] 

We  require  a  term  D  90(t)  which  is  not  the  result  of  visoous  friction. 
One  system  is  shown  in  Figure  7  employing  output  rate  feedback  by  means  of  a 
tachometer  generator  driven  from  the  output  shaft,  whose  output  voltage  is  pro¬ 
portional  to  output  speed. 


Thus:- 


J  aae0(t)  .  a  80(t) 

“  '  O  T  •  “ 


dt 


at 


e»(t)  =  e(t)  - 


dt 


dt 


Motor  torque  =  K  0'(t) 

a  e0(t) 


dt 


=  K  6(t) 


6o(t)  = 


_ K  6i(t) 

(j  T?  +  (f  +KH)  D  +  K 


The  effective  friction  coefficient  has  thus  been  increased  from  f  to 
(f  +  K  Kj ),  and  it  may  be  noted  that  the  additional  factor  is'  a  function  of 
quantities  capable  of  reasonable  regulation.  Other  more  subtle  means  are 
possible  for  obtaining  the  same  effect  such  as  extracting  the  back  e.m. f.  of 
the  drive  motor,  or  using  the  feedback  properties  of  a  falling  torque-speed 
characteristic  motor.  This  system  has  several  disadvantages,  tending  to 
increase  the  error  for  velocity  input  and  lowering  the  rate  of  response. 


4.6.2  Modification  to  the  controller  characteristic  - 
proportional  plus  error  derivative  control,  or 
"error  rate  damping",  and  addition  of  an  integral 
term 

The  preceding  discussion  has  emphasised  the  importance  of  the  error  - 
output  ohain  which  is  a  primary  factor  in  the  system  response  behaviour  and 
introduced  the  conception  of  a  controller  whose  performance  may  be  modified 
to  provide  the  desired  system  characteristic;  in  the  foregoing  the  controller 
may  be  taken  to  include  the  amplifier-motor  system  between  error  and  output. 

The  controller  forms  a  convenient  unit  in  which  to  adjust  the  system, 
particularly  if  it  includes  an  electronic  amplifier.  in  the  latter  case 
correcting  networks  may  be  embodied  which  can  be  readily  adjusted,  even  under 
operating  conditions  if  required. 

Our  previous  controller  characteristic  was:- 

Torque  =  K  0(t) 
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We  now  consider  a  generalised  characteristic  of  the  formy- 

.  .  d  6(t)  d89(t)  1 

Torque  =  j^K  9(t)  ±  l  — - ±  m  --^5-  J 


The  governing  equation  is  then: 


9(t)  .  _ <i Ht*), 

(J  Ds  ±  m  D8  +  fD  ±  ID  +  K) 

(J  D2  +  fD)  6j(t) _ 

[(J  ±  m)  D8  +  (f  ±  l)  D  +  K] 

Inspection  of  preceding  governing  equations  shows  that  f  must  be 
kept  to  a  minimum  to  reduce  power  wastage  and  the  ideal  arrangement  would  be 
to  have  f  =  0  and  l  set  positive  at  the  appropriate  value.  Reduction 
in  inertia  leads  to  improved  response  and  although  resulting  also  in  a 
reduction  in  5  can  be  compensated  for  by  adjustment  to  l.  Thus  m  must 
be  set  negative,  but  not  such  that  the  effective  inertia  is  reduced  to  a  very 
low  value,  when  instability  would  result. 

If  f  is  reduced  to  zero  the  steady  state  error  for  a  step  input  of 
velocity  u±  is  also  zero  and  the  transient  solution  is> 

e  ^Wnt  sin  _wn/l  -  £a.t'_ 

0  — r  =  —  . —  ■  g ssr  -  - 

U1  / 1  -  if 


6^  = 
ui 


which  compares  with  the  solution  for  viscous  friction  damping  alone 


where 


6(t)£  '  ^-/TT? 


f  sin  [%/ 


1  -  tf.t  +  <f> 


<p  =  tan 


-i  2g/l  -  S8 

2£8  -  1 


These  are  plotted  for  comparison  in  Figure  8.  Thus  the  addition 
of  derivative  of  error  known  also  as  "error  rate  damping"  provides,  a  ready 
method  of  adjusting  the  response  of  the  system,  although  the  second  deri¬ 
vative  of  error  is  difficult  to  obtain  and  seldom  used.  This  method  is 
free  from  the  disadvantages  mentioned  in  relation  to  output  rate  feedback. 

Hence  the  torque  equation  is  (in  the  absence  of  viscous  friction) 

Torque  =  K  9(t)  +  l 

i.e.  proportional.  +  derivative. 

The  action  may  be  visualised  by  consideration  of  the  action  of  the 
two  torque  components  after  application  of  a  step  input  of  velocity  u± 
to  &  system  with  X,  <  1 :  - 
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Immediately  after  application  of  the  input  step  the  error 
rate  is  equal  to  input  velocity  and  the  rate  term  provides 
a  maximum  assisting  torque. 

As  the  output  speeds  up  the  excess  of  input  speed  over  output 
speed  diminishes,  thereby  reducing  the  error  rate  torque. 

As  the  output  accelerates  further  it  overtakes  the  input  and 
the  error  rate  and  its  torque  become  negative,  thus  reducing 
the  overshoot. 

The  system  is  also  known  as  "anticipation  control",  but  this  is  really 
a  misnomer,  the  action  being  an  augmenting  or  boosting  action.  True  anti¬ 
cipation  requires,  for  example,  a  feed  forward  to  the  controller  from  the  input 
in  such  a  way  that  the  controller  can  take  action  before  the  generation  of  an 
error  or  derivative  of  error. 

The  system  is  additionally  known  as  a  "phase  advance"  system,  for  reasons 
which  will  be  considered  later. 

When  derivative  and  viscous  friction  damping  are  both  employed  the  res¬ 
ponse  characteristics  are  between  the  basic  responses  of  the  independent  systems. 

Presence  of  viscous  friction  re-introduces  the  problem  of  eliminating 
steady  state  errors  and  by  considering  preceding  sections  it  may  be  deduced 
that  these  may  be  eliminated  by  the  addition  of  an  integral  term,  thus  giving 
the  overall  controller  characteristic : - 

Torque  =  K  0(t)  +  l  ■  -  ^  --  +  K'^ 8(t)  dt 

T  I 

proportional  derivative  integral 

As  has  been  considered  previously  integral  control  tends  to  a  de¬ 
stabilising  effect  and  hence  the  various  parameters  must  be  carefully  selected 
for  the  best  solution. 

4. 7  The  stability  of  simple  systems  of  higher  order 

As  the  control  becomes  more  complex  it  is  not  easy  to  deduce  the  stability 
of  the  system;  the  systems  discussed  above  have  been  limited  to  seoond  order 
differential  equations  which  can  result  in  a  neutrally  stable  but  not  a  divergent 
oscillatory  response.  With  higher  order  systems  the  latter  is  a  possibility 
and  since  solution  of  the  governing  differential  equation  becomes  a  major  under¬ 
taking  there  is  clearly  need  for  a  short  cut  to  determination  of  the  system 
stability  which  does  not  require  a  complete  solution  of  the  differential  equa¬ 
tion. 


(a) 

(*) 

(c) 


It  may  be  deduced  from  previous  sections  that  a  system  may  be  described 
by  a  general  equation  of  the  form:  - 

(®n  D11  +  an-i  D""1  + .  +  ao)  «0(t)  -  **(*) 

and  that  the  stability  of  the  system  is  dependent  upon  the  transient  solution, 
i.e.  the  complementary  function.  The  stability  depends  then  on  the  roots  of 
the  auxiliary  equation  which  may  be  written:  - 

(an  vf1  +  an_i  irf1’1  + 


+  &o)  »  0 


I 


V 


r 

I 

?  ■ 

p :  * 


t 


I*- 

t 

■fe 
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or  (%  -  Yi )  (®a  -  Ya ) .  (%-*  -  Yn-i )  (%  ~  Yn)  3  0 

where  Yi »  Ya  .  Yk . Yn  8X6  the  roots  of  the 

auxiliary  equation. 

If  is  real  the  corresponding  solution  is  exponential;- 
0o(t)  =  A  e*kt 

and  is  damped  provided  Yic  is  negative. 

If  Yk  is  complex  the  roots  will  exist  in  pairs  of  the  form  (o^  t  j  Pk) 
giving  solutions  of  the  form:  - 

60(t)  =  A»eakt  00s (0kt  +  fe) 

which  w gain  is  damped  provided  is  negative. 

We  can  now  write  down  conditions  for  stability  based  on  the  coef¬ 
ficients  of  the  differential  or  auxiliary  equation:-/ 

(1)  The  real  parts  of  the  roots  of  the  auxiliary  equation 
must  be  negative. 

Thin  may  be  interpreted,  if  we  exclude  neutrally  stable  systems  such 
as  the  basic  second  order  control,  as;- 

(la)  All  powers  of  D  or  m  must  be  present. 

(lb)  The  coefficients  of  D  or  m  must  have  the  same  sign. 

(2)  The  coefficients  of  D  or  m  must  fulfil  certain  algebraic 

relationships . 

The  criterion  has  been  given  by  Routh  for  n  =  3:- 

8lj  fig  > 

and  when  n  =  4  there  must  be  satisfaction  additionally  of  the  relation:. 

ax  ag  a^  a^  a$  ^  +  as  ax 

It  is  important  at  this  juncture  to  appreciate  that  whilst  this 
criterion  establishes  whether  a  system  is  stable  or  not  it  gives  no 
indication  of  the  degree  of  stability. 

The  significance  of  the  above  is  indicated  by  consideration  of  the 
location  of  the  roots  in  the  m  plane,  when  it  will  be  seen  that  the  roots 
must  lie  in  the  left-hand  half  of  the  plane  for  stability,  that  no  conjugate 
roots  may  lie  on  the  imaginary  axis  and  that  no  multiple  order  sero  roots 
nmy  lie  at  the  origin  (Figure  9). 

Figure  9  also  serves  to  introduce  the  conception  that  there  are  only 
six  basic  forms  which  the  roots  of  the  characteristic  or  auxiliary  equation 
nmy  take,  of  which  only  two  are  basically  favourable  from  a  stability  view¬ 
point. 
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The  root  forms  are 
Root  1 
Root  pair  2 
Root  3 
Root  pair  4 
Root  5 
Root  pair  6 

It  should  be  emphasised  that  there  are  systems  with  multiple  loops  which 
may  have  roots  in  the  right-hand  half  of  the  m  plane  and  yet  are  stable,  pro¬ 
vided  the  main  loop  remains  closed.  The  standard  works  such  as  Reference  1 
should  be  consulted  for  further  consideration  of  these. 

4.8  Response  of  a  typical  system  to  a  harmonic  oscillation  of  input 

The  step  input  function  considered  so  far  has  represented  the  extreme 
and  generally  unrealisable  case  of  an  instantaneous  change  in  command  and 
although  the  resulting  response  of  output,  error,  etc.  has  been  considered  on 
a  time  basis,  no  consideration  has  been  given  to  the  response  of  the  system  on 
a  frequency  and  phase  basis  as  the  result  of  a  periodic  function  of  input. 

It  will  have  been  inferred  from  results  given  above  that  3ince  some  of 
the  systems  considered  have  exhibited  oscillatory  tendencies,  they  will  when 
subjected  to  an  appropriate  harmonic  input  exhibit  resonance  properties  leading 
to  various  amplitude  and  phase  shift  errors;  the  latter  may  of  course  result 
in  excessive  or  dangerous  errors  in  control  operation. 

4.8. 1  Periodic  and  harmonic  motion 

When  a  function  is  repeated  at  a  regular  time  interval  it  is  said  to  be 
periodic,  the  simplest  form  of  periodic  motion  being  a  sine  or  cosine  function, 
which  is  a  harmonic  periodic  motion  and  has  the  form:  - 

A  A 

0=0  sin  ut  or  0  =  0  cos  wt 

A 

Y/here  0  is  the  maximum  value  of  amplitude,  and  the  periodic  time  of  motion 


u 


where  T  is  the  reciprocal  of  the  frequency  of  oscillation  in  cycles  per  unit 
time. 


m 

= 

+ 

■a 

m 

= 

+ 

a  ±  jp 

m 

= 

0 

m 

= 

+ 

DP 

m 

= 

- 

a 

ra 

= 

- 

a  +  JP 

>  Unfavourable 


Favourable 


Now  the  significance  of  the  harmonic  motion  is  revealed  by  the  Fourier 
theorem  which  states  that  any  periodic  motion  of  circular  frequency  w  may 
be  represented  by  the  vector  summation  of  a  series  of  harmonic  functions  of 
frequency  u,  2 w,  3w,  etc,  and  appropriate  phase  relation,  thus  justifying 
the  use  of  a  sinusoidal  function  as  the  base  form  of  periodic  motion. 

Then: 

f(t)  =  Aq  +  AI  sin  (ut  +  fe)  +  Ao  sin  (2ut  +  + . 

where  \jr; ,  fg ,  etc.  are  the  appropriate  phase  angles. 
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The  Fourier  series  may  also  be  written: - 

f(t)  =  a!  sin  ut  +  aa  sin  2wt  +  ....  an  sin  nut  +  bo  +  bj,  oos  wt  + 
....  +  bn  cos  nut  +  ... 


where  Ana  =  an2  +  bns 


tan  i((n 


Many  highly  developed  techniques  exist  for  determination  of  the 
various  coefficients,  i.e.  for  performing  harmonic  analysis.  Mechanical 
harmonic  analysers  are  now  common  tools  in  industry. 

4. 8. 2  System  linearity  and  the  principle  of  linear 
superposition 

Much  servo-system  theory  is  based  upon  the  assumption  of  linearity 
in  the  system,  which  in  effect  means  that  the  input-output  relationship  is 
continuous  and  independent  of  the  magnitude  of  input.  Whilst  this  is  rarely 
realised  in  practice  it  serves  as  a  working  hypothesis  which  is  then  capable 
of  mathematical  analysis.  Thus  the  response  of  a  linear  system  to  a  har¬ 
monic  input  is  a  harmonic  output  of  identical  frequency,  but  differing 
possibly  in  amplitude  and  phase;  the  response  of  a  linear  system  to  a 
complex  input  is  thus  the  vector  sum  of  the  individual  responses  of  the 
system  to  the  harmonic  components  obtained  by  Fourier  analysis  of  the  input. 

This  is  then  an  example  of  the  principle  of  linear  superposition. 

The  Fourier  theorem  thus  provides  the  link  between  harmonic  and 
step  function  response. 

4. 8. 3  Zero  displacement  error  system  with  harmonic 
input 

For  the  zero  displacement  error  system  considered  above  the  govem- 

=  u n8  6i(t) 

fr>na  ei(t) 

(D3  +  2£  o>n  D  +  on*) 

=  B  cos  wt 

utjj8  B  cos  ut 
(2s  +  2?  oh  D  +  «n9) 


ing  equation  has  been  shown  to  be :  - 

(D2  +  2?  wn  D  +  wns)  0o(t) 

9o  (t) 

If  0i(t) 

0o(t) 


Now  the  solution  to  this  equation  will  take  the  form  of  a  transient 
solution  added  to  a  forced  solution  as  previously.  In  a  stable  system 
the  transient  solution  will  die  out  leaving  the  forced  solution,  i.e,  the 
particular  integral  only,  which  by  normal  operational  methods  gives:-. 
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60(t)  = 


/(i  -  aa  )a  +  (2£a)8 


cos  (ut  -  tan"1— ■? 


1  -  a8 


where 


This  solution  is  known  as  the  steady  state  response  to  a  harmonic  input 
and  must  not  be  confused  with  the  alternative  use  of  the  term  steady  state  to 
depict  a  condition  of  rest  or  equilibrium. 

The  above  solution  represents  an  output  vector  with  the  properties 


where 


Amplitude  =  —  — —  ,  .  _  -  3 

/(I  -  a2)2  +  (2£d)a 

Phase  s  -  tan-1  —  -2£d  ■  =  * 

(1  -  d2)  * 

©i(t)  =  B  cos  ut 

e0(t)  =  B  %  cos  (ut  +  ^r) 


Relative 
to  the 
input 
vector 


It  will  be  seen  that  in  the  region  d  =  1,  for  small  values  of  £,  Mj^  = 

Bg=a5r  ia  considerably  greater  than  unity,  hence  %  is  known 
/  (1  -  d8]8  +  (2£d)z  “ 

as  the  dynamic  magnifier.  As  £  -f  0  1%  ->  00  at  d  =  1,  and  *  =  --•  radians;  as 
d  -v  00  |  ->  -7t  radians.  ^  ' 

These  results  are  plotted  in  Figure  10  for  various  values  of  £  and  it 
will  be  noted  that  the  maximum  value  of  Mjj  occurs  at  a  value  of  d,  d  <  1. 

4.  8. 4  Zero  displacement  error  system  with  error  rate 


The  response  of  this  system  to  unit  harmonic  response  input  l.cos  wt  is 
given  by  the  particular  integral  of: 


Hence 


e0(t)  = 


Wn  cos  wt  -  2£w  wn  sin  wt 


(D2  +  2£  u^  +  una) 


(1  -  d8)2  +  (2£d)8 


•  sin  (ut  +  tan-1  2£d  -  tan-1  — ^ — 


1  -  d8 


where 


%  = 


1  +  (2£d)z 

(1  -  d8)3  +  (2£d)3 


i  =  -  tan-1  +  tan" 

1  -  d2 


Comparing  this  solution  with  that  for  the  system  without  error  rate 
damping  in  4.8.3,  it  will  be  noted  that  "the  phase  has  been  advanced"  by 

tan"1  2£d,  and  that  as  d  -+  co  -v  -  ~  . 
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These  results  are  plotted  in  Figure  1 1 . 

It  will  be  noted  that  Ui/L  =  1  for  all  values  of  g  at  d  =  /  2. 

5.0  Correlation  between  transient  response  and  harmonic  response 

It  will  be  noted  from  comparison  between  Figure  6,  and  Figure  10 
and  Figure  1 1 ,  that  there  is  a  close  relation  between  the  properties  of  the 
transient  response  and  the  harmonic  response,  and  considerable  significance 
attaches  to  this  in  the  design  and  improvement  of  control  systems.  For 
example  peaks  in  the  harmonic  response  of  a  complex  system  would  be  a  guide 
to  the  relative  damping  of  the  natural  modes  of  oscillation  of  the  transient 
response. 

As  a  general  value  My  =1.3  -  .1.4  is  often  considered  to  be  a 
reasonable  value  for  a  practical  system. 


5.1  The  use  of  complex  numbers  to  solve  for  functions 
of  time  and  frequency  - 


ejut  = 

cos  wt  +  j  sin  wt 

If  we  assume  0j,(t)  = 

l.sin  wt 

= 

Imaginary  j  e  j 

0o(t)  = 

A.  sin  (wt  +  ^r) 

= 

Imaginary  jAe^1*^ 

If  we  return  to  the 

governing  equation: 

(JD8  +  fD  +  K)  60(t)  =  K6i(t) 

By  rewriting  this  equation  in  complex  form:- 
(-Jwa  +  jwf  +  K)  Ae^ut  +  =  Kejut 

Aej* 


“n 


(-w3  +  2j  uZun  +  «na) 


Thus  we  may  write 


0 

where  ( j«)  is  a  vector  ratio  of  the  output  to  input  expressed  as  a 
function  of  the  complex  variable  jw. 

It  is  usually  more  convenient  to  work  in  non-dimensional  quantities, 

hence  writing  d  =  —  :  - 
«n 


AeJ*  = 


ef  <*•> 


“n 


[(j«)a  +  (jw)  +  una' 


6o 


1 


61  '  (1  -  d8)  +  2&)d 
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(1  -  d2) 


(1  -  a3)2  +  (2ga )* 


-  j 


2gd 


(1  -  a3)3  +  (2£L)a 


j  D 


This  represents  a  vector  of  magnitude  /C3  +  D3  having  a  phase  angle 


■  tan"4^-  (Figure  12). 


Thus:- 


Vector  magnitude 


/  (1  -  d2)3  +  (2gd)8 


Vector  phase  -  tan  1 


0O  .  . 

The  locus  of  the  vector  ■*—  ( jd)  is  plotted  in  Figure  12  for  various 

°i 

values  of  £. 

The  solution  as  a  function  of  time  is: 


fe(t) 


1 


sin  wt  -  tan 


1  -  d2 


i(1  -  d2)2  +  (2Sd)a 
i*ich  was  obtained  in  4.8.3. 

By  the  use  of  complex  numbers  we  have  been  able  to  obtain  a  vector  ratio 
of  output  to  input  as  a  function  of-  the  complex  frequency  variable  (ju)  or 
( jd)  which  describes  the  behaviour  of  the  system  when  it  is  subjected  to  a 
harmonic  input  in  the  absence  of  transient  terms;  i.e.  the  steady  state 
response  to  a  harmonic  input.  Under  these  conditions  if  we  write: - 


flo(t)  =  Q  (D)  6i(t) 

it  is  valid  to  substitude  ju  for  D  and  obtain  the  vector  ratio  of  the  complex 
variable :  - 


I”  (ju)  =  Q  (J«) 

It  is  a  simple  matter  to  obtain  from  this  relationship  the  solution  to 
the  original  governing  equation  as  a  time  function  in  response  to  the  harmonic 
input. 

6.0  The  conception  of  a  transfer  function 

Introduction  of  the  concept  of  a  transfer  function  has  been  deliberately 
avoided  in  the  earlier  part  of  the  discussion  to  prevent  confusion  resulting 
from  the  number  of  ways  in  which  it  may  be  defined;  at  the  commencement  of  the 
consideration  of  control  systems  the  system  appeared  as  a  number  of  elements 
which  combined  to  establish  a  differential  equation  from  which  the  response  was 
calculated.  Later  on  it  was  shown  that  control  systems  could  be  conveniently 
grouped  into  orders  and  the  conception  of  a  controller  was  introduced.  There 
was  thus  established  the  idea  of  an  input  to  an  element  which  was  operated  on 
within  the  element  to  provide  an  output  and  the  general  relation  obtained:- 
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e0(t)  =  q  (d)  ©i(t) 

where  Q  (d)  represents  an  operator  which  modifies  the  input  to  provide 
the  output.  This  quantity  is  called  the  transfer  function  of  the  element, 
(not,  it  my  be  noted,  the  transfer  function  of  the  system).  Alternatively 
we  may  transform  to  the  complex  form:- 


(jw)  =  Q  (jw) 

where  Q  ( jw)  represents  the  element  transfer  function  as  a  vector  based 
on  the  complex  variable  ju.  Reference  to  Equation  (7)  will  recall  that 
the  field  system  of  the  generator  and  amplifier  represented  a  simple  lag 

which  had  the  transfer  function  ( — — — j  and  that  in  general  the  transfer 

Vi  +  tD  1 

function  has  the  form  KS(D)  or  KG(  jw)  where  K  represents  the  gain,  or 
part  of  the  transfer  function  invariant  with  frequency,  and  G-(jw)  repre¬ 
sents  the  complex  part  varying  with  frequency. 

We  now  have  a  means  of  representing  the  behaviour  of  an  element  in  a 
concise  and  convenient  manner,  so  that  its  input-output  performance  is  always 
known;  this  is  particularly  convenient  in  a  complex  system  such  as  Figure  13, 
where 

e3 '  =  Qi.ei 

63  =  Qa*  63 

64  =  Q3.  63 

The  important  conclusion  which  follows  is  that  by  multiplying  transfer 
functions  we  can  obtain  the  overall  input -output  relation:- 


or 


04 


Qi  Qa  Qa  •  ®i 

K1G1  (  jw)  .K3G2  (  jw).KaG3  (jw) 
KG(ju) 


where  K  =  Ki.K9.K3 

©  (d«)  =  GiU«).^(dw).G»(j«) 

The  foregoing  concerns  elements  or  groups  of  open  chain  elements, 
but  our  main  concern  is  with  closed  loop  systems;  we  may  depict  this  in 
block  diagram  form  (Figure  13)  where  Q  =  Q1.Qj.Q3  etc.  depicts  the  element 
transfer  functions  grouped  together.  Now  if  a  signal,  say  e  ,  originates 
in  the  loop  it  will  be  transformed  on  arrival  at  the  same  point  after  tra¬ 
versing  the  system  as:- 


—  ^  Q1.Qj.Q3  etc.  |  c  —  —  Q  € 

where  Q  represents  the  loop  transfer  function  and  is  descriptive  of  the 
functioning  of  the  complete  system. 


So  that:- 
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0o(t) 

=  q(d)  e(t) 

and  since 

e(t) 

=  fli(t)  -  e0(t) 

e0‘(t) 

Q(D)  o-i 

i  +  q  (d)  01 

We  have  now  effected  a  subtle  transfer  of  emphasis  from  the  input-output 
behaviour  of  the  system  to  the  System  response  on  the  basis  of  a  loop  transfer 
function  which  relates  output  to  error  ( not  to  input). 


Thus:- 


r<*> 

e(jw) 


ef  <*•> 


KG(jw) 

6i(ju)  -  60(jio) 

k ajju) 

1  +  KG(ju) 


where  KG(  jw)  is  the  loop  transfer  function;  in  American  literature  this  is 
known  sinqply  as  the  transfer  function  of  the  closed  loop  system,  or  merely  the 
transfer  function. 


Since 


(j«) 


1 

1  +  KG(dti) 


we  obtain  the  important  characteristic  equation  for  the  closed  loop  system  with 
unity  feedback ;- 


1+  K3(jw)  =  0 


It  may  be  noted  in  passing  that  if  the  closed  loop  system  is  opened  at 
a  convenient  point  and  a  signal  sent  round  the  loop,  measurement  of  the  signal 
on  return  to  the  same  point  gives  a  measurement  of  the  loop  transfer  function. 
This  may  be  done  by  feeding  in  a  range  of  harmonic  signals  at  different  fre¬ 
quencies  and  measuring  the  dynamic  magnifier  and  phase  angle  of  the  returning 
signal  relative  to  the  input ^  thus  enabling  the  harmonic  response  locus  to  be 

plotted,  i.e.  the  locus  of  —  (ju);  the  latter  function  is  called  the  harmonic 

response  function. 


7.0  Use  of  the  Laplace  transformation 


The  Laplace  transform  method  is  an  operational  calculus  approach  which 
enables  impulsive  and  discontinuous  functions  to  be  dealt  with,  provides  the 
transient  and  steady  state  solutions  simultaneously  and  is  readily  able  to  take 
account  of  initial  conditions.  The  method  is  superior  to  the  classical  operator 
technique  for  the  solution  of  high  order  equations,  and  by  use  of  the  transform, 
conversion  of  the  differential  equation  in  the  time  domain  to  a  function  of  a 
complex  variable  can  be  effected  in  which  this  function  can  be  manipulated  by 
the  normal  laws  of  algebra;  this  functional  relationship  can  be  returned  to 
the  time  domain  by  inverse  transformation. 
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The  method-  is  used  in  American  literature  almost  exclusively,  but 
al though  it  is  a  powerful  method  of  considerable  value,  its  use  in  ele¬ 
mentary  problems  may  be  no  advantage  over  the  classical  methods  already 
discussed.  Due  to  substitutions  vdiich  can  be  effected  when  certain 
restrictions  are  applied  to  the  system,  casual  inspection  of  American 
literature  may  lead  to  considerable  confusion  in  the  mind  of  the  reader  as 
to  the  use  of  the  method  and  even  in  the  definition  of  properties  of  the 
system  such  as  transfer  functions.  Indiscriminate  substitution  of  the 
Laplace  parameter  may  also  lead  to  error  unless  initial  conditions  are 
reconciled  and  a  true  appreciation  of  the  method  obtained.  It  is  for 
this  reason  that  consideration  of  the  method  has  been  left  to  this  stage 
in  the  paper. 

A  function  of  time  and  of  the  complex  variable  are  related  by  defini¬ 
tion  as:- 

f(p)  =  t  f(t) 

Where  f(p)  is  a  function  of  the  Laplace  parameter  p* 

p  is  a  complex  variable  of  the  form  a  +  jw 
i l  denotes  application  of  the  transform  integral 
f(t)  is  a  function  of  time. 

The  transform  integral  is:- 

00 

f  e”?*  f(t)  dt  =  £ 

' o 


There  are  certain  limitations  to  unrestricted  application  of  the 
above  which  do  not  concern  us  here. 

If  we  consider  a  unit  step  function  of  time;- 

f(t)  =  1 


f(p) 


Hence 


f(p) 


■i 

—  for  unit  step  function. 
P 


For  an  exponential  decay: - 


f(t)  =  e”wt 


*  The  Laplace  parameter  p  is  denoted  by  s  in  American  literature 
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t  (e-ut)  =  F e"(P  +  w)t  dt  =  — 

Jo  p  +  w 


For  a  sinusoidal  function:- 


£  (sin  <Jt)  =  ^  sin  wt  e dt 

oo 

=  SJ  { (eJut  - e'pt  dt 

l  (sin  ut)  =  jr  (  — “ : - " — ~  )  =  —  - — a 

\  p  -  JW  p  +  0“  /  p  +  u 

Some  frequently  used  transform  pairs  are  given  in  Table  III, 

It  follows  that:- 

£^K-f(t)J  =  K  f(t)  =  K  f (p) 

£  )  +  fa('t)J  =  fi(p)  +  ^(p) 

£  [  ft  f(t)j  =  P  f(p)  -  f(0  +) 

t  £_  f(t)l  =  pa  f(p)  -  pf(0  +)-|  (o  +) 

L  dt  J 

-  ir  *p[/f(t)  at]t .  (o  +) 

The  terms  f  (0  +)  may  occasion  some  surprise  but  is  necessary  due  to 
the  use  of  discontinuous  functions  and  means  the  value  of  f(t)  as  t  approaches 
zero  from  the  positive  direction.  Consider  for  example  a  unit  step  function 
occurring  at  t  =  0  then  at  t  =  (0  -),f(0  -)  is  zero,  but  at  t  =  (0  +), 
f(0  +)  is  unity. 

Various  theorems  and  applications  of  the  initial  and  final  value  theorems 
enable  the  relationships  to  be  manipulated  and  short  cuts  effected  (Reference  2). 

Inverse  transformation  denoted  by  ~1  enables  functions  of  p  to  be 
converted  back  to  the  time  domain  e.g. 


( j? .  «■ ) " 


sin  ut 
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We  return  now  to  consideration  of  the  simple  system  denoted  by 
Equation  (11 )  in  the  absence  of  load  torque:- 

D  (JD  +  f )  60(t)  =  K  0(t) 


Since 


%  ■/} 


22  %  =  j 

D  (D  +  25  Uft)  90(t)  =  a*2  6(t) 

V*  e0(t)  +  25  ft*  D  e0(t)  =  Un2  6(t) 

We  can  now  take  the  Laplace  transform  of  each  term  to  obtain: - 

1  ]  =  ^  6°^  "  p  9°  (°  +)  “  (°  +) 

£  25  o)n  ~°-J  =  25  wa'-p  0o(p)  -  0O  (0  +) 

t  J^2  e(t)J  =  Una  0(p) 

Thus  60(p)  = - r~* - ~ - -1  0(p) 

p  p  +  25  wnJ. 

- ‘  V' 

transfer  function 

jp  6p  (0  +)  +  (0  +)  +  25  Wn  6o  (0  +)J 

P  [p  +  25  WnJ 


initial  condition  transform 


If  we  consider  60(p)  in  the  form:- 


0  (p)  =  we  may  eventually  resolve  the 

B(p) 


expression  by  partial  fractions  into  the  form:- 


B(p)  (P  -  Yi )  (p  -  Ya ) 


+  • .  • . 


a(p)  a(p) 

provided  - -  is  proper  and  contains  no  repeated  roots,  or  if  ■  is 

B(p)  B(p) 

improper  or  contains  repeated  roots  then  we  proceed  by  the  normal  rules  of 
partial  fractions.  The  inverse  transform  may  then  be  obtained  for  each 
fraction  individually,  considerable  lose  being  made  of  the  property:- 


parameter,-  or:- 
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Since  the  Laplace  transform  method  is  favoured  in  the  U.S.A.  the 
transfer  function  is  defined  there  as:- 


KGf  )  -  transform  of  the  output 

'  '  ~  Laplace  transform  of  the  error 


which  if  encountered  without  the  background  discussed  above,  puts  an  appar¬ 
ently  complicated  construction  on  a  basically  simple  function. 


Analysis  and  synthesis 


Opportunity  will  now  be  taken  to  make  an  appreciation  of  the  material 
covered  so  far  in  relation  to  control  system  performance. 


In  the  first  part  of  the  discussion  the  system  existed  as  an  assembly 
of  rather  unrelated  elements  and  solution  of  the  governing  differential 
equation  yielded  the  transient  response  as  a  time  function  relation  between 
output  and  input.  The  conception  of  order  of  control  was  introduced  and 
it  was  established  that  the  majority  of  systems  could  not  have  the  perfect 
response  but  were  a  compromise.  The  presence  of  friction  and  inertia  in 
a  practical  system  was  shown  to  raise  the  order  of  the  differential  equation 
and  to  introduce  oscillatory  behaviour  depending  upon  the  value  of  damping 
present.  Consideration  of  a  typical  system  employing  viscous  friction 
damping  alone  was  shown  to  be  exceedingly  inefficient,  and  alternative 
methods,  of  introducing  damping  were  considered  3uch  as  output  rate  feedback, 
and  a  derivative  of  error  term  or  error  rate  damping;  the  addition  of  a 
negative  second  derivative  term  to  reduce  the  effective  inertia  wa3  also 
considered  together  with  an  integral  term  for  removal  of  steady  state  error. 
This  introduced  the  conception  of  a  controller  between  error  input  and  out¬ 
put,  comprising  for  example  an  amplifier  and  motor  to  provide  the  required 
circuits  and  power  amplification  for  obtaining  the  desired  response  of  out¬ 
put  to  error.  Consideration  was  given  to  determining  the  stability  of 
simple  systems  without  the  necessity  for  a  complete  solution,  by  inspection 
of  the  coefficients  of  the  governing  equation  and  the  important  point  was 
made  that  Routh's  criterion  did  not  determine  the  degree  of  stability. 

System  response  to  a  harmonic  input  was  examined  and  the  correlation  exist¬ 
ing  between  the  form  of  response  to  this  and  to  step  functions  of  input 
noted.  The  block  interpretation  of  systems  led  to  the  idea  of  an  element 
transfer  function  defining  completely  the  input-output  function  of  the 
element  and  introduced  the  important  loop  transfer  function  of  a  closed 
loop  system  relating  output  to  error,  which  contained  all  the  information 
required  about  the  system.  Finally  the  vector  locus  and  Laplace  transform 

were  introduced  and  the  identity  of  ^  ,  D,  jo,  and  p  established  for  a 

system  subject  to  a  harmonic  input  in  which  the  initial  conditions  were 
zero;  the  characteristic  equation  of  the  closed  loop  system  with  unity 
feedback  ms  shown  to  be  1  +  KGr(p)  =  0.  It  was  also  shown  that  the  loop 
transfer  function  oould  be  established  by  opening  the  loop  at  any  point, 
injecting  a  signal  and  noting  the  modifications  to  it  introduced  in  its 
passage  back  to  the  starting  point;  if  the  signal  comprised  a  harmonic 
input  whose  frequency  could  be  varied,  it  was  possible  to  determine  the 
loop  response  by  measurement  of  the  dynamic  magnification  and  phase  of  the 
returning  signal;  from  this  the  harmonic  response  locus  could  be  plotted. 


The  crux  of  the  control  system  problem  has  been  shown  to  lie  then 
not  in  an  analysis  of  the  performance  of  a  typical  system  in  the  form  of 
the  input-output  response,  but  in  the  combined  problem  of  establishing  the 
design  of  a  system  which  will  have  the  required  response  and  then  ensuring 
that  its  response  will  in  fact  be  the  one  desired.  This  is  the  problem  of 
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synthesis  and  much  of  the  confusion  which  a  newcomer  to  the  subject  experiences 
is  his  inclination  to  work  on  the  basis  of  analysis,  when  the  wealth  of  control 
system  technique  is  devoted  to  synthesis; '  this  results  in  his  being  confronted 
with  apparently  abstract  quantities  such  as  the  vector  locus  of  the  complex 

variable  jw ,•— ( jw) ,  rather  than  input-output  functions  such  as  0o(t)  and 

8j(t)*  It  may  appear  as  if  solution  of  the  differential  equations  or  algebraic 
treatment  of  the  roots  of  1  +  KG(p)  =  0  must  enable  synthesis  to  be  effected, 
but  in  the  case  of  a  complex  system  this  involves  vast  amounts  of  computing, 
and  although  the  effect  of  changes  in  the  system  can  be  calculated,  it  is  virt¬ 
ually  impossible  to  obtain  the  optimum  arrangement  of  the  complete  control 
system  except  by  chance.  Considerable  effort  has  been  expended  on  analysis  of 
the  transient  response,  which  for  many  years  was  the  only  tool  available  and  a 
number  of  techniques  and  graphical  methods  developed  for  easing  the  determination 
of  system  response  and  stability. 


Modem  methods  of  synthesis  are  founded  primarily  on  a  frequency  response 
basis,  using  the  correlation  between  transient  and  harmonic  response  to  ensure 
that  the  transient  response  will  be  satisfactory.  Much  of  the  technique  is 
based  upon  properties  of  the  loop  transfer  function  and  its  vector  locus  £G(jw); 
an  alternative  version  of  the  latter  known  as  the  inverse  transfer  function 
1—1 

locus  — t — r  or  KG  (jw)  is  sometimes  more  convenient.  A  major  part  of  the 
KG(j<d) 

work  comprises  adjusting  the  K  function  and  modifying  the  G(  ju)  locus  until 
the  desired  performance  is  achieved.  Hie  degree  of  stability  is  also  established 
by  these  techniques. 


One  major  advantage  of  this  approach  is  that  by  dealing  in  vector  loci, 
system  elements  may  be  dealt  with  on  a  transfer  function  basis  by  determining 
the  element  transfer  function  experimentally,  even  though  the  differential 
equation  of  it  cannot  be  formulated  mathematically.  This  method  lends  itself 
readily  to  the  use  of  analogues;  if  for  example  the  complex  machine  or  engine 
is  to  have  an  electrical  controller,  the  transfer  function  may  be  determined 
experimentally  and  an  electrical  analogue  designed  to  give  the  same  response 
as  the  actual  engine.  Calculations  and  synthesis  can  proceed  on  the  basis  of 
the  measured  response  and  checks  made  on  the  actual  control  apparatus  employing 
the  analogue  rather  than  the  actual  machine. 


9.0  The  design  process 

The  first  requirement  in  the  design  is  a  specification  defining  the 
response  under  certain  basic  conditions  which  include :- 

(1)  The  degree  of  stability  required  following  a  step  or  pilse 
disturbance. 

(2)  The  error  0(t)  which  is  acceptable  for  specified  variations 
of  input  0j_(t)  or  load  Ti/t). 

(3)  The  phase  shift  between  input  and  output  which  is  allowable 
and  the  frequency  band  width  to  be  covered. 

With  a  knowledge  of  the  particular  function  of  the  control  system  the 
designer  may  then  be  able  to  lay  down  the  basic  properties,  e.g.  order  of 
control  etc. 


The  procedure  then  comprises: - 


(4)  Construction  of  a  block  diagram  representing  the  physical 
elements. 
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(5)  Formulation  of  the  element  transfer  functions  by 
calculation  or  experiment;  where  non-linearities 
exist  the  system  must  be  linearised  by  the  accepted 
mathematical  methods. 

(6)  The  system  may  then  be  analysed  and  the  response 
determined,  thus  completing  the  analysis. 

(7)  Should  there  be  undesirable  features  present  in  the 
response,  the  techniques  of  synthesis  must  be  applied 
and  the  system  adjusted. to  give  the  required  response. 

It  is  not  the  purpose  of  this  Report  to  consider  the  design  of 
control  systems  except  in  so  far  as  this  helps  to  clarify  the  engine 
response  fate  problem,  but  some  consideration  must  be  given  to  the  basic 
properties  of  the  transfer  function  locus  on  the  complex  plane,  i.e. 
KG(ju)  if  the  picture  is  to  be  complete. 

10,0  Some  basic  properties  of  the  KG(.jm)  locus  and  its  inverse 

locus  KG-M  j&i) 

It  has  been  established  that  for  the  closed  loop  system  with  unity 
feedback  defined  by:- 


6i(jw)  -  0o(jw)  =  e(jw) 

KGr(ju) 


e, 

0i 


TT<*) 


I -<d“) 


1  +  KG(ju) 

KG(ju) 

1  +  KG(jw) 

KG^Cjw) 

1  +  KG(  jai) 

KG(ju) 

=  1  ♦  IffiPQw) 

e 

10.1  Graphical  determination  of  jw)  from  the  KG(jw)  locus 

A  representative  locus  of  KG( jw)  is  given  in  Figure  14»  the  vector 
being  of  amplitude  OB  and  phase  Point  A,  which  will  be  shown  to 

be  very  important,  is  at  -1  on  the  real  axis,  i.e.  the  (-1  +  Oj)  point. 

Then  AB  represents  the  magnitude  of  the  vector  1  +  KG(  jw),  and 
phase  -e. 


te  .  «■»(*.)  .  Vju) 

AB  1  +  KG(jo)  ' 


Thus 
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and  the  phase  angle  =  angle  GBA  =  <j>  -  e. 

The  locus  of  3—  (jw)  may  -then  be  re-plotted  (Figure  14)  and  it  will  be 
°i 

recalled  that  the  magnitude  of  this  vector  is  the  dynamic  magnifier  %.  It 
may  be  noted  that  the  triangle  ABO  represents  vector  AO  s  0,  vector  OB  a  0o 
vector  AB  =  0j_  to  an  appropriate  scale. 

10.2  Stability  and  the  (-1  +  Oj)  point,  the  Nyquist  diagram 

When  the  locus  passes  through  the  (-1  +  Oj)  point  it  follows  that  a  signal 
originating  in  the  loop  will  pass  round  the  system  and  appear  at  the  starting 
point  as  a  signal  of  equal  amplitude  and  in  phase  with  it.  The  signal  is  then 
self-sustaining  without  the  need  for  an  input.  Thus  the  (-1  +  Oj)  point  corres¬ 
ponds  to  the  case  of  poles  lying  in  the  imaginary  axis  when  discussing  roots  of 
the  characteristic  equation  plotted  in  the  m  or  p  plane;  a  locus  enclosing 
the  (-1.+  Oj)  point  and  representing  a  divergent  oscillation  corresponds  to 
roots  of  the  characteristic  equation  lying  in  the  right  hand  half  of  the  m  or 
p  plane. 

The  criterion  for  stability  may  thus  be  stated: - 

If  the  KG(jw)  locus  is  traversed  from  the  u  =  0  point  in  the  w  =  00 
direction,  the  system  is  stable  provided  the  (-1  +  Oj)  point  is  always  to  the 
left  of  the  locus. 


This  definition  is  valid  for  single  loop  systems  and  for  multi -loop  systems 
with  stable  inner  loops. 

The  criterion  is  named  after  Nyquist  who  first  stated  it  in  relation  to 
amplifier  design  and  it  is  of  course  well-known  in  relation  to  the  performance 
of  feedback  amplifiers. 

The  harmonic  response  locus  KG(jw)  is  often  known  as  a  Nyquist  diagram, 
but  it  must  be  observed  that  the  latter  term  is  subject  to  frequent  abuse  and 
is  often  applied  indiscriminately  to  polar  response  plots  based  on  functions  of 
jw  other  than  the  KG(ju)  locus. 

Further  consideration  of  the  relationship  between  the  three  basic  loci 
is  given  in  Appendix  I. 

10.3  KG(j<j)  and  the  roots  of  the  characteristic  equation 


It  is  convenient  at  this  juncture  to  refer  back  to  the  transient  response 
via  the  Laplace  parameter  and  functions  plotted  on  the  complex  plane  and  to 
consider  again  the  roots  of  the  characteristic  equation.  The  major  difficulty 
in  dealing  with  the  techniques  used  in  synthesis  is  to  obtain  an  understanding 
of  the  interrelation  between  different  approaches  to  a  given  problem  and  to  put 
a  practical  interpretation  on  some  of  the  more  abstract  aspects.  This  is 
especially  true  when  dealing  with  the  KG(jw)  locus  for  the  first  time. 


It  has  been  shown  that  p 
a  complex  number. 

Thus  we  may  write:- 

.  '  ( 

P  = 
|  f(p)  = 

i.e.  f  (a  +  ju)  = 


is  a  complex  number  and  that  f (p)  is  also 


a  +  ju 
u  +  jv 
u  +  jv 
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arid  the  value  of  u  +  jv  varies  continuously  with  changes  in  a  and  w. 

Thus  f(p)  may  be  plotted  (Reference  1 )  in  terms- of  a  and  u  to 
form  a  network  of  constant  a  and  constant  w  lines  (Figure  15). 


Since 

f (a  +  jw) 

= 

u  +  jv 

and 

P 

= 

a  +  jw 

22 

3a 

= 

1  “*!£  = 

j 

3f 

3p 

= 

df 

dp 

3f 

3a 

= 

1! 

• 

af 

dp 

= 

f'(p)  . 

3f 

3W 

= 

2£  .  22  = 
3p  3u 

■3  • 

df 

dp 

=  jf'(p) 

Thus 

the  spacings 

of 

the  constant  a 

and 

w  contours  are  equal  at  any 

point  and  intersect  orthogonally. 

Reference  to  Figure  15  shows  also  the  plot  of  constant  a  and  to  lines 
on  the  p  plane  which  was  considered  in  Figure  9.  Again  these  are  equally 

spaced  and  cut  orthogonally,  hence  the  f  (p)  plane  is  a  conformal  trans¬ 
formation  of  the  p  plane. 

We  have  already  seen  that  the  system  response  is  governed  by  the 
characteristic  equation: - 

1  +  KGr(p)  =  0 

or  1  +  f(p)  =  0. 

The  solution  to  this  equation  is  the  value  of  p  which  gives 
f(p)  =  -1.  This  may  be  determined  directly  in  the  form  (a  +  jw)  by  read¬ 
ing  from  Figure  15  the  values  of  a  and  to  at  which  the  network  inter¬ 
sects  the  (-1  +  Oj)  point.  The  graphical  work  may  be  reduced  by  plotting 
the  value  of  f  (a  +  jw)  for  a  =  0,  and  constructing  the  complete  net  by 
the  orthogonal  properties  of  the  a,  w  contours. 

It  follows  that  f(a  +  jw)  for  a  =  0  is  f(jw)  and  that  the  line 
for  a  =  0  corresponds  to  the  harmonic  response  locus,  when  f(jw)  =  K&(jw). 
If  thi3  locus  passes  through  the  (-1  +  Oj)  point  it  clearly  corresponds  to 
a  root  of  the  characteristic  equation  with  a  zero  real  part  and  when  it 
passes  to  the  left  of  the  (-1  +  Oj)  point  it  corresponds  to  a  positive 
real  part  to  the  complex  roots  discussed  previously,  and  henoe  to  a  diver¬ 
gent  oscillation.  This  provides  additional  clarification  of  the  signi¬ 
ficance  of  KG(jw). 

It  is  possible  to  extend  the  graphical  treatment  given  above  to 
obtain  a  complete  solution  to  the  transient  response  to  step  inputs 
(Reference  1). 
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10,4  Constant  dynamic  magnifier  and 
~  (jw)  in  the  KG(jw)  plane 


ise  contours  of 


Wq  /  \ 

Contours  of  constant  1%  and  phase  for  “  (ju)  may  be  plotted  in  the 
KG  plane  (Figure  16). 

« _  _ _ 


Since  %  is  the  magnitude  of 


1  +  KG(ju) 


this  may  be  re-arranged  with  x  and  y  as  Cartesian  co-ordinates  of  the  locus 
of  KG(  ou)  to  give 


y.  ,  =  . . 

L  ^  ■ 1 J  ■  (%I  ■ 1)a 

Thus  the  locus  of  constant  Mjj  is  a  circle  whose  centre  is  located  on  the 
negative  real  axis  with  co-ordinates:- 


V  Hm  - 1  / 


and  has  a  radius :- 


It  will  be  noted  that  this  locus  is  independent  of  u,  thus  providing  a 
generalised  plot. 

Similarly  when  \|r  is  the  phase  lead  of  output  on  input  we  may  write:- 

N  =  tan  ljr 

l"W)'  ♦(r-aj  =  °-25  p-jfr] 


The  centre  of  the  circle  is  at  co-ordinates:- 

(-0-5'  *s) 

and  has  radius :- 

If  we  consider  a  plot  of  the  K&(ju)  locus  with  %  and  N  circles 
superimposed  it  is  clear  that  we  can  deduce  the  system  response  from  it  readily 
and  if  required  plot  conventional  curves  of  and  ^  against  frequency 
(figure  17).  The  three  curves  show  the  effect  of  an  increase  in  K:  - 

(a)  Larger  My  values 

(b)  Less  phase  shift 

(c)  No  ohange  in  basic  locus  shape 
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Further  increase  in  gain  would  result  in  a  violently  oscillatory 

system. 

It  is  seen  that  adjustment  of  gain  may  enable  a  particular  1%  or 
t  specification  to  be  met,  but  that  gain  adjustment  alone  cannot  always 
meet  both  requirements. 

10.5  Adjusting  K  in  the  G(.ju)  plane  for  a  selected  value 


If  we  draw  the  locus  of  G(jw)  i.e.  K  =  1  we  oan  draw  a  line  at 

4 

an  angle  6  to  the  real  axis  (Figure  18),  such  that  6  =  sin”4  rr-  .  A 
circle  can  then  be  dram  so  that  the  sin"1  rr-  line  and  the  locus  Gr(  jw) 
are  tangents.  This  is  possible  since:-  .  . . 


sin  6 


where  ^  is  the  scale  factor  of  the  drawing  in  so  fair  as  the  true  value 
of  %  is  concerned. 


Now 


OC 

OB 


~  where  BC 


is  a  perpendicular  to  QAfram  B 


OC 


1 

K 


This  gives  the  required  value  of  K  for  the  criterion 

seleoted  and  the  locus  of  KG(jw)  oan  then  be  drawn  as  described  prev¬ 
iously.  The  resonant  frequency  oan  be  determined  as  the  value  of  u  at 
which  the  locus  is  a  tangent  to  the  M)[  circle. 


This  section  provides  a  further  illustration  of  the  value  of  work¬ 
ing  in  terms  of  the  K£(ju)  locus  rather  than  the  overall  response 

(ju)  since 


1  +  KG(jw) 
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0 

Any  change  in  K  would  necessitate  a  complete  replot  of  the  ( jw) 
locus,  where  a  a.  it  merely  requires  a  change  of  scale  for  the  K&(  jw)  locus. 

10.6  Gain  and  -phase  margin 


These  terms  are  often  used  in  design  work  and  may  be  stated  in  slightly 
different  ways;  one  definition  is:- 

Gain  margin  The  reciprocal  of  the  KG(  jw)  magnitude  for  a  phase  of 
-x  radians.  Typical  value  2,5  to  10. 

Riase  margin  The  angle  between  the  negative  real  axis  and  the  KG(jw) 
vector  ^or  the  frequency  at  which  the  magnitude  of  KG(  jw)  is  unity.  Typical 
value  30°  to  60° . 

10. 7  The  inverse  locus  KG-1 ( jw) 

Similar  properties  exist  for  the  locus  of  KG-1 ( jw)  as  for  the  locus 
of  KG(jw)  discussed  above. 


Since 


- -  =  1  +  KG_1  ( jw) 

eT  (j“») 


a  very  simple  vector  relation  exists  on  the  KG-1  plane  (Figure  18)  and 

constant  %  contours  are  circles  centred  at  (-1  +  Oj)  and  radius  —  ;  constant 

“M 

N  contours  are  diameters  of  the  %  circles  making  an  angle  tan  N  with  the 
real  axis,  where  N  =  tan  f  and  V  is  the  phase  lead  of  input  on  output. 


10.8 


Adjusting  K  in  the  G-1  ( jw)  plane  for  a  selected  value 
of 


The  process  is  similar  to  that  used  on  the  G(jw)  plane.  The  sirT1  tr* 

1  » 

line  is  drawn  for  the  desired  where  sin  o  =  —  and  a  circle  drawn  such 

that  the  locus  and  the  sin"'1  line  are  tangents  to  it.  Since  the  centre  of 

the  circle  is  located  at  (-1  +  Oj)  in  the  KG  plane  it  follows  that  if  the  locus 
is  plotted  to  its  own  scale,  then  OA  =  K  the  desired  value  of  .K, 


10.9  Modifying  the  G(jw)  locus 


It  has  been  shown  that  modification  to  the  gain  K  alone  cannot  usually 
produce  the  required  performance  and  the  G(  jw)  function  itself  must  be  re¬ 
shaped.  This  may  be  effected  by  modification  to  the  system  arrangement,  or  by 
introduction  of  "shaping  circuits"  which  are  in  effect  active  or  passive  elements 
designed  so  that  the  introduction  of  their  transfer  functions  into  the  loop 
produces  phase  and  amplitude  changes  in  the  response  locus.  The  addition  of 
phase  advance  for  example  to  a  lagging  system  whose  K&(  jw)  locus  passes  through 
the  (-1  +  Oj)  point  can  move  the  locus  into  the  stable  region.  . 

The  general  subject  of  reshaping  the  G(  jw)  locus  is  however  outside 
the  scope  of  this  paper. 
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10.10  The  relative  merits  of  the  KG(ju')  and 
loSI 

It  may  be  difficult  to  appreciate  the  relative  merits  of  the  two 
methods  of  plotting  from  the  above,  since  the  properties  discussed  appear 
to  be  very  similar  for  either  locus.  Since  the  low  frequency  end  of  the 
locus  lies  near  the  origin  for  the  inverse  locus,  this  method  of  plotting 
is  preferable  for  consideration  of  the  low  frequency  response,  but  where 
the  high  frequency  response  is  concerned  this  is  best  dealt  with  by  the 
Nyquist  diagram  or  KG(jw)  locus.  The  methods  are  thus  complementary. 

The  inverse  locus  however  possesses  a  particular  advantage  in  deal¬ 
ing  with  multiple  loop  feedback  systems,  or  systems  containing  dynamic 
elements  in  the  feedback  path.  This  may  be  illustrated  by  reference  to 
the  output  rate  feedback  system  considered  in  4.6. 1.  It  was  shown  that 
the  basic  equation  was  of  the  form:- 

6'(t)  =  Q(t)  -  Q(D)60(t) 


In  the  absence  of  the  additional  feedback  term  we  may  write 


©o(t)  =  KGi{D)0(t) 


Hence  with  the  feedback  term:- 


or 


0o(t) 

e0(t) 

—  M 


KGi (D)  0(t)  -  Q(D)0o(t) 

_JSd£L_  .  e(t, 

1  +  Q(D)  KGi (D) 

KSx(dM) 

1  +  Q(ju)KGj.(jw) 


Thus  with  the  harmonic  response  function  the  locus  must  be  re¬ 
calculated  for  any  change  in  the  feedback  term  Q( jw). 

The  inverse  locus  however  is;- 


f  U)  =  KGi^Ou)  +  Q(ju) 

which  enables  the  additional  loops  to  be  included  by  simple  vector  addition 
of  the  feedback  term  to  the  inverse  locus  of  the  system  which  is  to  be 
modified. 

11.0  The  logarithmic  or  asymptotic  form  of  locus 

The  asymptotic  form  of  locus  was  proposed  over  twenty  years  ago  and 
has  been  intensively  developed  in  the  field  of  oonmunioations  engineering 
during  recent  years. 

The  technique  comprises  plotting  the  magnitude  and  phase  of  the 
transfer  fuxvotion  separately'  against  the  frequency  in  logarithmic  co¬ 
ordinates:  some  surpris-  may  be  occasioned  by  this  since  it  appears  to  be 
a  reversion  to  techniques  considered  earlier  and  rejected  as  unsuitable  for 
synthesis.  The  significance  lies  however  in  the  use  of  logarithmic  co¬ 
ordinates. 
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The  logarithm  of  the  magnitude  of  KG(  ju)  at  any  point  may  be  written: 

log  KG(ju)  =  log  K  +  log  G( ju) 

The  phase  or  angle  of  KG(ju)  is:- 
Ang  KG(ju)  =  Ang  G(ju) 


For  a  complex  system; - 


KG(ju) 


KiG^ju)  •  KQ»(ju) 


At  a  particular  value  of  KG(ju)  we  may  write  for  the  amplitude 
log  KG(ju)  =  log  Kx  +  log  Ka  -  log  Ka 

+  log  Gi(ju)  +  log  Gg ( ju) 

-  log  G3 ( ju) 

/ 

Thus  former  multiplication  and  division  become  addition  and  subtraction 
and  change  in  gain  merely  requires  the  addition  or  subtraction  of  a  constant. 

Since  this  approach  is  allied  to  general  techniques  of  a  similar  nature 
in  communications  engineering,  the  same  units  are  employed  i.e.  decibels.  These 
originated  largely  as  the  result  of  the  human  ear  having  a  basically  logarithmic 
response.  Thus  successive  increases  in  sound  intensity  of  say  ten  times  would 
create  the  impression  of  equal  increments  to  the  human  ear,  rather  than  equal 
multiplication. 

It  follows  that  if  a  power  level  of  P3  is  established,  greater  than  Px  f 

the  increase  in  level  of  Pa  over  Px  is  log  =—  relative  to  a  datum  at  Px .  Since 

hi 

it  is  conventional  to  work  in  log10  we  can  write 

Pa 

Change  in  power  level  =  logj.0  in  bels 

hi 


e  10  logjo  in  decibels, 
hi 

Now  if  the  input  and  output  impedances  are  equal  and  say  R: 


Hence  change  in  voltage  or 

signal  amplitude  level  =  20  lo&0 


Where  V0  and  Vi  are  the 
output  and  input  voltages  or 
amplitudes 


decibels. 


When  V_  =  2  Vi  the  change  in  level  is  20  .  0.  .501  a  6  decibels  (db)  and 
it  will  be  noted  that  a  further  doubling  of  the  output  would  result  in  a  further 
increase  of  6  db  making  12  db  in  all  from  the  datum  level. 
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It  is  not  uncommon  for  the  subject  to  be  discussed  on  a  similar 
basis  in  respect  of  frequency,  in  that  octaves  are  used  to  describe  change 
in  frequency  level  or  pitch;  thus  doubling  the  frequency  represents  an 
increase  in  frequency  level  of  one  octave. 


Hence  change  in  frequency 
level  relative  to  a  datum  at 

It  is  convenient  to  abbreviate  the  expression: - 

20  log!  o  KG(ju)  =  20  logio  K  +  20  logxo  G(ju) 

to  Lm  KG(ju)  where  the  units  are  decibels.  The  phase  relation  Ang  G(ju) 
is  usually  in  degrees. 

As  an  example,  we  may  now  consider  the  transfer  function:- 

■  TT* 

Lm  K&( ju)  =  Lm  K  +  Lm  (1  +  ju)  1 


=  Lm  (1  +  ju)  1 

The  amplitude  and  phase  logarithmic  plots  are  given  in  Figure  19. 

The  amplitude  plot  exhibits  interesting  properties  which  must  be  considered. 


When  u  -*  oo 

Lm  G(  ju)  -+•  -  oo 

u  -♦  0 

Lm  G(ju)  -»  0 

u  =  0.5 

I*m  G(ju)  =  -  1  db 

u  =  1.0 

Lm  G(ju)  =  -  3  db 

u  =  2.0 

Lm  G(ju)  =  -  7  db 

u»>1 

Lm  G( ju)  o  log  (ju)  1 

Thus  at  the  higher  frequencies  the  amplitude  falls  by  6  db  for  each 
doubling  in  frequency;  this  is  often  quoted  as  the  amplitude  falling  6  db 

per  octave. 

Thus  two  asymptotes  describe  the  function,  one  at  a  constant  level 
of  0  db  and  the  other  starting  at  0  db,  u  =  1  and.  falling  at  6  db  per 
octave;  if  desired  the  three  corrective  points  calculated  above  may  be 
used  to  remove  the  error  in  the  region  u  =  1. 

The  point  of  intersection  of  the  two  asymptotes  is  known  as  the 
break  point,  which  here  is  at  u  =  1;  this  is  known  as  the  corner  frequency 
and  occurs  at  the  value  of  u  which  makes  the  real  and  imaginary  parts  of 
the  function  numerically  equal. 


/fa  \ 

=  logg  I  -s-  )  octaveB  where  fa  and 
'  1 '  f !  are  the  respective 
frequencies. 
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It  will  be  noted  that  variation  in  gain  factor  K  merely  changes  the 
reference  level  of  the  whole  plot. 

The  other  basic  factors  of  the  transfer  function  can  be  treated  in  a 

similar  manner  and  the  general  method  forms  a  powerful  tool  in  synthesis,  comple¬ 

mentary  to  the  vector  loci  techniques. 

12.0  General  consideration  of  transfer  function  loci  shapes 

We  have  seen  that  the  characteristic  equation  may  be  factorised  to  give;- 

(P  -Yi  )  (P  -  Ya  )  •  •  •  (P  "  Yk)  •  •  •  (P  “  Yn-i )  (P  -  Yn)  =  0 

and  that  the  roots  maybe  real,  imaginary,  or  in  complex  pairs,  and  in  general 
will  have  the  basic  forms 

(P  +  0) 

(PT  +  1 ) 

(P2  +'2£  Wn  p  +  Uj!8) 

By  making  a  substitution  of.  the  form  u  =  f(p)  we  obtain  basic  forms:- 

ju 

(ju  +  1) 

(ju)a  +  2?  (ju)  +  1 

and  in  general  the  transfer  function  will  consist  of  factors  of  the  form;- 


(ju)11 
(ju  +  I)11 

.  [(ju)9  +  2£  (ju)  +  1 

(juX  +  Ip 

[(juX)3  +  25  (juX)  +  1  11 

Where  KG(ju)  ha3  the  form  K  •  etc. 

C(ju)  •  D( ju) 


One  of  the  most  important  properties  of  the  transfer  function  locus  lies 
in  the  fact  that  an  individual  skilled  in  the  art  can  by  visual  inspection  alone, 
obtain  a  considerable  amount  of  data  from  it  relating  to  the  general  properties 
of  the  system  or  element. 

In  Figure  20  a  collection  of  the  more  common  and  useful  loci  are  tabled, 
together  with  a  diagram  of  the  locus  shape  for  each.  Provided  the  transfer 
function  locus  is  available  for  a  system  or  element  the  latter  may  generally  be 


-  49  - 


Report  No.  R.187 


classified  even  if  the  exact  mathematical  function  is  unknown.  Reference 
4  may  be  consulted  for  a  comprehensive  consideration  of  harmonic  loci. 


Amongst  the  functions  tabled  are  the  basic  first  and  second  order  of 
control  functions  in  (l)  and  (2),  the  ideal  derivative  term  controller 
function  in  (6)  and  the  ideal  integral  terms  controller  function  in  (3). 
Since  the  latter  require  infinite  gain  at  infinite  and  aero  frequency  res¬ 
pectively,  they  are  impracticable  for  general  use  and  would  in  fact  be 
undesirable  in  most  systems  due  to  the  presence  of  noise  and  random  signals 
which  require  the  operating  frequency  range  to  be  limited  to  the  minimum 
acceptable  band  width.  The  functions  in  (7)  and  (4)  are  more  represent¬ 
ative  of  realistio  derivative  and  integral  term  responses. 


Typical  shaping  circuits  for  obtaining  various  input-output  element 
characteristics  are  shown  in  Figure  21. 


An  important  operator  considered  in  4.5.1  was  the  finite  delay  and 

11 

it  was  shown  to  have  the  transfer  function  form  -swr  or  — rr  .  The  de- 


«TD 


stabilising  effect  was  discussed  and  may  be  illustrated  by  consideration  of 
the  effect  of  this  operator  on  the  harmonic  response  locus. 

Now  =  cos  u  -  j  sin  u 


and  represents  a  rotation  of  the  input  vector  by  an  angle  -  u,  3ince  the 
amplitude  is  unchanged.  Thus  at  any  point  on  the  harmonic  response  locus 
frequency  ujj,  the  delay  operator  rotates  the  vector  and  oauses  it  to  lag 
by  u  angular  units  (Figure  22). 

13.0  Conclusions 

The  properties  of  elementary  control  systems  have  been  considered, 
the  basio  terminology  described  and  the  conception  of  order  of  control' 
introduced. 

The  presence  and  effect  of  lag  terms  in  a  practical  system  was  noted, 
particularly  that  of  a  finite  time  delay  and  the  conception  of  a  controller 
whose  response  could  be  adjusted  to  give  the  desired  response  was  introduced. 

The  importance  of  the  properties  of  the  roots  of  the  characteristic 
equation  were  emphasised  and  the  location  of  them  in  the  m  plane  noted  in 
relation  to  stability  and  system  response.  The  damping  term  in  a  simple 
system  was  shown  to  result  in  three  basic  forms  of  the  roots  of  the  charac¬ 
teristic  equation,  leading  to  underdamped,  critically  damped  and  overdamped 
responses. 

Initially  step  functions  of  input  were  employed,  but  consideration  of 
a  harmonic  input  illustrated  the  amplitude  and  phase  properties  of  a  system 
with  variation  in  input  frequency  and  the  limitation  of  response  with  fre¬ 
quency  dictated  by  the  constants  of  the  system.  The  correlation  between 
transient  and  harmonic  response  was  emphasised  and  the  ability  to  transform 
from  one  domain  to  the  other  by  various  techniques  noted.  It  was  shown 
that  consideration  of  the  harmonic  response  after  the  decay  of  the  initial 
condition  terms,  reduced  the  response  function  to  the  particular  Integral 
solution  of  the  differential  equation  and  use  of  complex  numbers  indicated 
that  the  oomplex  frequency  parameter  ( jw)  could  be  substituted  for  D  or 
m  to  give  the  harmonic  response  on  a  vector  basis. 
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The  conception  of  an  element  and  loop  transfer  function  was  introduced 
relating  output  to  input  for  the  element  and  output  to  error  for  the  loop,  which 
comprised  complex  numbers  completely  describing  the  properties  of  the  element  or  ' 
system  under  dynamic  conditions.  A  brief  note  was  recorded  on  the  use  of  the 
Laplace  transformation  method,  which  is  much  favoured  in  the  U.  S.A.  and  is  a 
powerful  method  for  the  solution  of  complex  systems. 

The  crux  of  the  problem  was  reached  in  considering  analysis  and  synthesis, 
the  transient  approaoh  leading  to  analysis  of  system  behaviour  but  being  relatively 
intractable  for  determining  the  optimum  system  arrangement  to  meet  the  desired 
performance  specification  i.e.  the  process  of  synthesis, 

Achievement  of  a  satisfactory  input-output  response  is  obviously  the  goal 
of  the  designer,  but  it  was  shown  that  achievement  of  the  desired  response  by 
synthesis  was  powerfully  aided  by  the  frequency  response  properties  of  the  output- 
error  function  to  a  harmonic  input  i.e.  the  loop  transfer  function  KG  expressed 
as  a  function  of  (ju)  or  of  a  non-dimensional  frequency  parameter  ( ju)  and 
known  as  a  Nyquist  diagram.  Various  graphical  constructions  were  considered 
leading  to  establishment  of  the  stability  and  gain  performance  with  frequency 
and  of  the  relation  between  plots  in  the  KS(  jw)  plane  and  the  roots  of  the 
characteristic  equation.  The  relative  merits’  of  the  harmonic  response  locus 
KG(jw)  and  the  inverse  locus  KG~1(jw)  were  considered,  the  latter  being  of 
advantage  when  dealing  with  multiple  loops  and  dynamic  elements  in  the  feedback 
path.  The  logarithmic  locus  was  mentioned  briefly  as  a  further  aid  to  synthesis 
and  finally  a  number  of  the  more  commonly  used  elements  and  their  transfer  function 
loci  were  tabled.  In  this  connection  the  de-stabilising  effect  of  a  finite  time 
delay  was  illustrated  by  demonstrating  its  effect  as  an  operator  on  the  harmonic 
response  locus. 

Consideration  of  the  engine  response  rate  programme  may  raise  doubts  as 
to  the  necessity  for  such  an  extensive  consideration  of  control  system  theory, 
particularly  techniques  of  synthesis  which  lie  rather  outside  the  scope  of  the 
programme.  Much  of  the  currently  available  literature  however  originated  in 
the  U.S.A.  where  the  fullest  exploitation  of  servo-system  technology  in  the  gas 
turbine  response  rate  field  has  been  attempted.  Whilst  this  is  unquestionably 
a  valid  approach,  it  has  resulted  in  some  of  the  basic  papers  which  concern 
fundamentally  simple  systems  being  stated  in  terms  which  to  a  non-specialist 
render  them  incomprehensible ,  but  which  would  be  readily  understood  if  expressed 
in  a  more  conventional  manner.  For  this  reason  Part  I  of  the  paper  has  been 
extended  to  cover  most  of  the  terminology  used  in  current  papers. 

It  may  be  concluded  from  consideration  of  the  wealth  of  servo-system 
technique  already  available,  that  provided  the  engine  element  functions  can  be 
defined  in  an  acceptable  manner  and  the  system  lags  and  non-linearities  esta¬ 
blished,  solution  of  the  engine  problem  by  established  means  should  be  possible 
(Reference  3). 

Where  however  unpredictable  changes  such  as  variable  lags  or  general 
inconsistency  are  present,  or  the  system  cannot  be  linearised  by  acceptable 
means,  it  is  inevitable  that  the  problem  will  become  exceedingly  complex. 
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AEPEMDIX  I 


Further  consideration'  of  the  relationship  between 

6 

the  loci  KG(jw),  K G_a(ju)  and  —  (jw) 

Further  consideration  of  the  basic  vector  diagram  may  help  to  clarify 
the  relationship  between  the  three  primary  loci  K&( jw),  KG’1  ( jw)  and  ^2  (jw) 
i.e.  the  Nyquist  diagram,  the  inverse  locus  and  the  harmonic  loou a. 

D.  B,  Welboum  has  drawn  the  attention  of  the  author  to  an  illustra¬ 
tion  based  on  a  particular  example.  Consider  for  example  a  damped  simple 
springe-mass  system: - 

Og D*  +  TXD'+  1)  60(t)  =  0i(t) 
where  T*  =  2£  Ta 

IShen  the  system  is  subjected  toharmonic  input s- 
(!§!)“  +  TjD  +  1)  0o(t)  =  A  cos  wt 

the  solution  for  the  output  will  have  the  form:- 

0o(t)  s  B  cos  (ut  -  f  ) 

By  writing  D  a  jw 

[t§  ( j«)a  +Ti(jw)  +  l]e0(t)  =  6i(t)  =  A  cos  wt 

It  is  now  possible  to  plot  a  vector  diagram  representing  this  equation, 
noting  that  the  operator  j  rotates  the  vector  through  9<f  and  that  j 
rotates  the  vector  through  18C f.  This  is  given  in  Figure  23(a).  UV,  Vtf 

and  IX  represent  the  magnitudes  of  0O,  wTi0o  and  -w*T§0o  respectively. 

These  vectors  sum  to  give  TJX  denoting  the  vector  d±.  It  follows  that 
VX  represents  the  difference  vector  0*  -  0O  i.e.  the  error  0,  If  the 
diagram  is  scaled  to  make  UV  a  1,  VX  represents  the  vector  defining  the 
0 

locus  (jw)  i.e.  the  inverse  locus. 

If  the  diagram  is  now  rearranged  (Figure  23  (b))  with  the  origin  at 
U  and  TJX  =  1,  vector  UV  defines  the  locus  (jw),  or  in  brief  the 
harmonic  locus. 

Finally  by  taking  the  origin  at  V  and  making  VX  *  1 ,  vector  VU 
0 

defines  the  locus  ^  (jw),  i.e.  the  Nyquist  diagram. 
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APPENDIX  II 


a  a  constant 

=  a  constant  coefficient  in  a  differential- equation 
=  a  Fourier  coefficient  e.g.  A^8  =  an8  +  bn2 
=  a  constant 
=  a  Fourier  coefficient 

=  the  real  part  of  a  complex  root  pair  (a  ±  Jb) 

=  a  constant 
=  a  constant 
=  a  Fourier  coefficient 

=  the  imaginary  part  of  a  complex  root  pair  (a  ±  Jb) 

=  a  constant 

s  the  real  part  of  (0  +  JD) 
a  a  constant  of  integration 
a  the  differential  operator  e  ^ 

=  the  imaginary  part  of  (C  +  JD) 

=  the  dimensionless  frequency  ratio  ~ 

=  the  natural  base  of  logarithms  =  2.7183 
=  instantaneous  e.m.  f. 

*  the  coefficient  of  viscous  friction 
s  a  function  e.g.  of  time  f(t) 

=  a  periodic  frequency  in  cycles  per  unit  time 

=  the  coefficient  of  viscous  friction  for  critical  damping  *  2t/jK 
=  the  first  derivative  of  a  function  e.g.  ^ 

=  the  value  of  a  function  as  the  .zero  is  approached  from  the 
positive  direction 

=  the  frequency  dependent  part  of  a  transfer  function  KG(  ) 

=  the  acceleration  due  to  gravity 
=  instantaneous  electrical  current 
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APPENDIX  II  (cont'd) 

polar  moment  of  inertia 
/“ 

a  constant 

a  stiffness  in  load  per  deflection  units 

a  gain  factor  e,g.  the  frequency  invariant  part  of  a  transfer 
function  K&(  ) 

un 

the  velocity  constant  -gg 
a  constant 

a  stiffness  in  load  per  deflection  units 
electrical  inductance 


the  coefficient  of  the  first  derivative  term  in  a  generalised 

(cLO  cl^  0  \ 

IC  +  &  +  m  J 

dynamic  magnifier  i.e.  ratio  of  dynamic  amplitude  to  equivalent 
static  deflection 


a  root  of  the  auxiliary  equation 

the  coefficient  of  the  second  derivative  term  in  a  generalised 

,  , ,  Z"  ( v  ,  a©  ■  da e  \ 

torque  equation  T  =  ^  K  +  £  +  m  ^3-  j 

the  phase  parameter  tan  \[r 

electrical  power 

the  Laplacian  parameter 

the  Laplace  transform  of  f(t)  i.e.  l  |f (t)J 
a  function  of  an  operator  or  complex  variable  e.g.  Q(d), 

Q(j«) 

electrical  resistance 

the  Laplacian  parameter  in  American  literature 
a  finite  time  delay  ' 
the  periodic  time  — 


instantaneous  time 


a  dimensionless  frequency  parameter 

the  real  part  of  the  complex  number  representing  f(p),  i.e. 
(u  +  jv) 

electrical  voltage 

the  imaginary  part  of  the  complex  number  representing  f(p), 
i.e.  (u  +  jvj 
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KG-(D) 

KG(p) 

KG(ju) 

KG(ju) 

y 

KG-1  (  )  = 


P 

r 

r 

6 

r  _ 

e  I- 

s . 


e 

n 

-c 

* 


w  = 

&>n  = 

Un  /  1-$f  = 

l  ~ 


APPENDIX  II  (cont'd) 


a  transfer  function 


the  inverse  form  of  a  transfer  function 
the  real  part  of  a  complex  root  e.g.  (a]c  +  j  Pk) 
the  real  part  of  the  complex  number  p,  i.e,  (a  +  jw) 
the  imaginary  part  of  a  complex  root  e.g.  (ak  +  j  Pk) 
torque 

a  factored  root  of  the  auxiliary  equation 


a  signal  originating  at  a  point  in  a  closed  loop  system 
a  phase  angle 

the  ratio  of  actual  damping  coefficient  to  the  coefficient 

for  critical  damping  =  ■s—  =  — =» 

2  /JK 

angular  displacement 
the  gain  of  an  amplifier 
a  time  constant 
a  phase  angle 

a  phase  angle  e.g.  in  a  Fourier  series  =  tan' 


■i  *2. 


frequency 

the  natural  frequency  of  ■undamped  oscillation  = 


the  natural  frequency  of  damped  oscillation 
denotes  application  of  the  Laplace  transform 


c 


Suffices 

=  critical  in  relation  to  the  damping  coefficient  f 
=  electrical  field 
=  input 

=  input  at  steady  state 
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APPENDIX  II  (cont'd) 

k  =  the  k'kh  term 

L  =  load 

o  =  output 

oss  =  output  at  steady  state 

M  =  motor 

n  =  the  n*'*1  term 

n  =  natural  frequency  in  relation  to  frequency  w 
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